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Schematic of the simple "shape-subtraction" 
paradigm for correction of myopia. Ri and R? are initial 
are final radii of curvature; to is the may i mum depth of the 
laser cut; S is the diameter of the optical zone. Reprinted 
from Munnerlyn et al 

(57) Abstract: Biomcchanical responses of the eye are used to improve photorefractive procedures. LASIK or PRK treatments, for 
example, can be improved by taking prc-opcrative measurement and predicting the corneas' s biodynamic response to the ablative 
treatment Predicitve use is mad of the biodynamic response of the cornea due to laser or mechanical keratectomy, that is, creating a 
corneal flap characterisitc of LASIK. Comparison of pre-flap and post-flap, (pre-ablation) data of the cornea such as corneal thick- 
ness, flap thickness, corneal topography and wavefront, for example, can provide predictive information applicable to modifying 
an ablation alorithm before the laser is engaged, either for a current operation or the development of a model. Modeling, by finite 
element analysis or other mathematical techniques, can also be used to predict post-operative outcomes based on a pre-operative (no 
flap cut or other surgical intervention) data for the cornea that is input for an accurate eye model that, in consideration of biodynamic 
response via the model, provides predictive information for optimizing the sucess of the refractive surgery and ultimately patient 
satisfaction. 
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METHODS AND INSTRUMENTS FOR IMPROVING 



REFRACTIVE OPHTHALMIC SURGERY 



Field of the Invention 



The invention relates generally to the field of vision correction and more particularly 
5 to methods and apparatus involving biodynamics^iomechanics of the eye to ultimately 
provide improved vision to patients through refractive correction. 



The future of laser refractive surgery is tremendously exciting as we move into the 
new millennium with the potential for ever-improved post-operative visual performance. In 

10 the past, the operative goal has been 20/20 visual acuity with zero residual refractive error. 
Criteria for a "successful" procedure were operationally broadened to 20/40 or better with ±1 
diopter of residual refractive error. However, neither surgeons nor future patients will likely 
be satisfied with such gross data of visual performance. Hence, considerable research effort 
is being devoted to develop "customized" procedures for each individual patient. The new 

15 ultimate operative goal is 20/10 visual acuity with aberration-free post-operative vision. Yet, 
how can these lofty goals be accomplished? First, lasers have been improved with the 
development of scanning, small spot systems, as opposed to the first generation broad-beam 
systems. The scanning systems have brought customized procedures into the realm of 
feasibility. Second, more comprehensive and sophisticated input data can be used to "guide" 

20 the laser based on individual patient data, as opposed to the simple sphere and cylinder of the 
past. Two types of approaches are currently being pursued - wavefront-guided and 
topography-guided procedures. Early results are promising, yet neither approach has 
demonstrated consistently superior results to non-guided procedures in controlled, scientific 
studies. Is a piece of the perfect vision puzzle still missing? Is there an additional, 

25 complimentary approach to "customization" that has yet to emerge? Are we leaping into the 
future ahead of our understanding of how these procedures may be implemented to be 
optimally successful? 

Approximately 60 percent of Americans have refractive errors, and millions of people 
are myopic worldwide. Many thousands of laser refractive surgeries are performed every 
30 year for the correction of myopia. These procedures will ultimately affect a huge number of 
individuals around the globe, and yet the corneal response to laser ablation is not well 
understood. Of the many individuals treated, about 15-50% of them do NOT achieve 20/20 



Background of the Invention 
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vision, which translates into very large numbers when the extreme popularity of refractive 
surgery is considered. It is crucial that the number of patients who achieve their targeted 
vision be increased in order to improve the overall quality of vision in this vast group of 
people. 

5 In the quest for the ideal corneal ablation, a major obstacle yet to be overcome is the 

inadequacy of current models for predicting the corneal response to ablation, which impacts 
the success of both topography-guided and wavefront-guided customized procedures. In its 
current state, the process of ablation design relies on a "black box" approach, wherein 
enormous efforts are directed at linking "input" variables (the ablation algorithm) to "output" 

10 variables (refractive error, visual acuity, glare, aberration, patient satisfaction, etc.), while the 
physical interface between input and output — the actual mechanisms of interaction between 
all relevant components — is largely neglected. Ablation and surgical outcome are ultimately 
linked on two levels: a deterministic, cause-effect relationship dictated by physical reality, 
and a statistical relationship, rooted in probability and defined in retrospective regression 

1 5 analyses of the same variables in large-scale clinical trials. The ability to objectively measure 
performance of the eye, as opposed to simple corneal shape, is critical to "customizing" 
ablation algorithms and generating an overall improvement of visual outcomes after 
refractive surgery. Anterior corneal surface topography cannot take into account 
contributions of optically important structures inside the eye, such as the posterior corneal 

20 surface and the crystalline lens. If a laser were programmed strictly with anterior topography 
data, the correction would be at best incomplete, and at worst simply wrong. Therefore, 
wavefront analysis is clearly important, . particularly if the ultimate goal is to correct higher 
order aberrations along with the sphere and cylinder. However, the question must be asked 
whether wavefront analysis alone is sufficient to fully predict visual outcomes. Will it 

25 replace corneal topography in the quest for the perfect "aberration-free"guided procedure? 
Or, on the other hand, is there a piece of the puzzle still missing? If so, can corneal 
topography complement the wavefront data to help complete the picture of corneal response? 

Initial attempts at photorefractive keratectomy (PRK) used a simple model presented 
by in 1988 by Munnerlyn et aL (Munnerlyn CR, Koons SJ, Marshall J. Photorefractive 
30 keratectomy: a technique for laser refractive surgery, J Cataract Refract Surg. 1988;14:46- 
52) The cornea was modeled as two refracting surfaces with a bulk material in between with 
a known index of refraction. In the treatment of myopia, the goal was to increase the anterior 
radius of curvature, thus decreasing the curvature of the anterior surface, as illustrated in 
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Figure 1. A simple geometric formula resulted, which assumed the targeted corneal shape 
was a function only of the ablation profile. This can be thought of as the "shape-subtraction" 
paradigm, based on a geometric approach to tissue removal and secondary curvature change, 
where the final corneal shape is thought to be determined solely by how much tissue is 
5 'subtracted' by the laser. Essentially, the cornea is a piece of plastic to be sculpted into the 
"ideal" surface shape. 

The equations described by Munnerlyn et al still serve as a starting point for 
developing current ablation algorithms. For a myopic ablation, the pre-operative cornea is 
modeled as a sphere of greater curvature than the desired post-operative cornea, which is also 

10 modeled as a sphere. The apex of the desired post-operative cornea is displaced from the pre- 
operative cornea by the maximal ablation depth, which is determined by the ablation zone 
size. The intervening tissue is simply removed or "subtracted" to produce the final result. 
This is illustrated in Figure 2 for a myopic, as well as a hyperopic profile which is similar in 
concept, but with a different planned ablation profile that produces increased post-operative 

15 curvature. This concept will be referred to as the "shape-subtraction" model of refractive 
surgery and treats the cornea as if it were a homogeneous structure, like a piece of plastic. 
The incomplete shape-subtraction paradigm permeates current thinking in refractive surgery, 
and forms the basis for both topography-guided and wavefront-guided procedures. 

The Munnerlyn approach, combined with empirical experience, has been relatively 
20 successful in correcting spherical and cylindrical errors for the majority of the patient 
population to date. However, a substantial number of individuals do not achieve 20/20 vision 
(only 50 - 85% do achieve 20/20, depending on who is reporting), despite greater than 90% 
being satisfied. In addition, significant post-operative optical aberrations are produced using 
conventional algorithms, hence the move toward developing aberration-reducing ablation 
25 profiles. 

After acceptance of PRK, laser in-situ keratomileusis (LASIK) was proposed as an 
alternative procedure. LASIK involves cutting a thin corneal flap of tissue, ablating within 
the stromal bed, and retaining the flap. LASIK has become increasingly popular since it is 
associated with less pain and faster visual recovery. In addition, it was thought that the 
30 preservation of the corneal epithelium in LASIK would minimize regression. Ultimately, 
however, both procedures produce similar visual outcomes. 
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When these procedures failed to produce the expected refractive outcomes with 
consistency, it became clear that the model was incomplete. Growing appreciation for the 
complexity of the corneal response coupled with a rapidly growing body of clinical 
experience led to a more empirical approach to algorithm design. Current ablation algorithms 
5 are proprietary to the particular manufacturer, and it is not known how similar they are to the 
original formula. However, it is believed that the ablation profiles have been empirically 
altered based on experience with real clinical data. 

Among a large array of possible input variables, the ablation algorithm is the most 
exquisitely controllable. Consequently, it is often the sole focus of refractive control in 

10 developing sophisticated corneal ablation routines based on topographic and/or wavefront 
analysis. However, the model of corneal response that contributed to unpredictability in 
noncustomized procedures still influences the refractive outcome and overall visual 
performance. Although empirically modified ablation algorithms have produced reasonable 
results in many patients, they are unlikely to provide the individual predictability that 

15 physicians desire and patients increasingly demand. Moreover, algorithms derived from 
these probabilistic models were optimized for the mean population response rather than that 
of the individual, and a certain degree of prediction error is inescapable. The current 
practical approach to excimer laser refractive surgery has become one of iteration between 
ablation recipe and surgical response, with most surgical decisions being made at the 

20 empirical level and confounded by an incomplete working knowledge of the cornea's 
physical behavior during and after ablation. Individual surgeons use their own "fudge" 
factors, making systematic interpretation of results difficult. Thus, the fundamental challenge 
to the future of custom corneal ablation is to develop deterministic models that can be 
successfully applied to the individual patient, rather than relying heavily on empirical 

25 modification. Important biomechanical determinants of corneal response to excimer laser 
ablation and a rationale for inclusion in customized corneal ablation models, presented 
herein, was previously unrecognized in ablative procedures. 



The cornea has been modeled as a biomechanical structure for refractive procedures 
30 that do not rely on excimer laser ablation, such as radial keratotomy. (Roy et al., 1996; Hanna 
et al., 1992;Vito et al. 1989; Bryant et al., 2000; Pinsky and Datye, 1991a,b) In these existing 
numerical models for refractive surgery, the cornea is assumed to be a solid material and the 
effect of the surgically-induced structural change on corneal deformation, produced by an 
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incisional or other procedure, is investigated. The only mechanical loading condition used in 
these simulations is the intraocular pressure applied to the posterior surface of the cornea. 
Living human cornea, however, is highly porous and filled with a biological fluid. In fact, 80 
% by weight is considered to be water. Also, the corneal endothelium serves as a biological 
5 "fluid-pump", maintaining a state of relative dehydration with a negative intra-stromal 
pressure (approximately 50 -60 mmHg below the atmospheric pressure) called the 
"imbibition pressure". In addition, the microstructure of the corneal stroma consists of 
approximately 300 ~ 500 lamellar layers throughout its thickness. Each layer consists of 
collagen fibers supported by a mstrix consisting of "proteoglycans" and water. When living 

10 human cornea is incised or ablated, the long collagen fibers are severed in a series of layers. 
In addition to the loss of the tensile force in the remaining collagen segments, the water in the 
neighborhood of the damaged portion of the cornea can no longer maintain the same level of 
negative ("imbibition") pressure. This change in water pressure inside the damaged portion 
of the cornea plays a crucial role in corneal deformation. However, this factor has never been 

1 5 considered in any of the existing simulation methodologies for refractive surgeries. 

The biomechanical response - as used herein, the term 'Tnomechanical response," 
sometimes used interchangeably with, "biodynamic response," means a mechanical or 
physical response to a perturbation or other stimulus - has not been directly incorporated into 
current ablation algorithms. As a result, predictability of the visual result currently 

20 achievable, even in "customized" wavefront-guided or topography-guided procedures, has 
been inadequate. There are significant differences in the material properties of living human 
corneas based on age, sex, race, climate, years of contact lens wear, etc. Thus, there are 
significant differences in the response of various individuals to similar ablative profiles. The 
ability to predict the biomechanical response of the cornea to a photoablative procedure, and 

25 account for it in ablation algorithm design would significant enhance quality of vision 
produced. Knowledge of individual biomechanical properties, prediction of individual 
reponses, and individual adjustment of generalized ablation profiles would allow individual 
customization, rather than optimizing to the mean population response. Therefore, a 
biomechanical model with predictive capabilities would significantly enhance the overall 

30 success rates of these refractive surgeries. 

Summary of the Invention 

The invention pertains most generally to improving visual quality through refractive 
surgery. An aspect of this general goal is to obtain the highest level of subjective satisfaction 

5 
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of visual quality as measured objectively. The gist of the invention resides in the predictive 
capability of the biodynamics of the eye. A wide variety of photoablative procedures can 
benefit from this predictive ability, in a number of ways. For example, LASIK or PRK 
treatments can be improved by taking pre-operative measurements and predicting the 
5 corneas's biodynamic response to the ablative treatment. In another embodiment of the 
invention, predictive use is made of the biodynamic response of the cornea due to laser or 
mechanical keratectomy, that is, creating a corneal flap characteristic of LASIK. Comparison 
of pre-flap and post-flap, (pre-ablation) data of the cornea such as corneal thickness, flap 
thickness, corneal topography and wavefront, for example, can provide predictive 

10 information applicable to modifying an ablation algorithm before the laser is engaged, either 
for a current operation or the development of a model. Modeling, by finite element analysis 
or other mathematical techniques, can also be used to predict post-operative outcomes based 
on a pre-operative (no flap cut or other surgical intervention) data for the cornea that is input 
for an accurate eye model that, in consideration of biodynamic response via the model, 

15 provides predictive information for optimizing the success of the refractive surgery and 
ultimately patient satisfaction. LASIK, PRK and other photoablative techniques benefit from 
these embodiments of the invention. 

In accordance with the foregoing, an embodiment of the invention is directed to a 
method for performing refractive corneal surgery. The steps include obtaining a pre- 

20 operative (e.g.., pre-flap cut for LASIK) diagnostic measurement data of the eye, determining 
ablation specifications based upon the pre-operative diagnostic measurement data, perturbing 
the cornea of the eye such that there is a biodynamical response, obtaining a post-pertubation 
diagnostic measurement of the eye, determining an adjusted ablation specification based at 
least in part on a comparison of the pre-operative and post-pertubation data wherein these 

25 data are indicative of a biodynamical response of the eye, and ablating a portion of the 
cornea using the adjusted ablation specification. As will be evident to one skilled in the art, 
there will be a link between the measurement comparison and the biodynamic response. For 
example, the post-pertubation measurement may indicate a greater or lesser degree of central 
corneal flattening and greater or lesser degree of peripheral corneal steepening and 

30 thickening. These structural and physiological changes are manifestations of the magnitude 
of the biodynamic response of the cornea to the pertubation. The magnitude of the 
biodynamic response can then be used to adjust the ablation specifications, e.g., less central 
pulses and more peripheral pulses. 
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Another embodiment is directed to a method for establishing a photoablative 
specification by obtaining a pre-operative (e.g. pre-flap cut) diagnostic measurement data of 
the cornea, detennining ablation specifications based upon the pre-operative diagnostic 
measurement data, perturbing the cornea of the eye such that there is a biodynamic response, 
5 obtaining a post-pertubation diagnostic measurement of the cornea, and determining an 
adjusted ablation specification based at least in part on a comparison of the pre-operative and 
post-pertubation data wherein these data are indicative of a biodynamical response of the 
eye. 

Associated with the immediately foregoing embodiment is a method for correcting 
10 and, preferably, optimizing existing algorithms for photoablative ophthalmic surgery that 
includes taking data of a statistically significant number of corneas. Biomechanical response 
has not been directly factored into current ablation algorithms. With a proper biodynamic 
model, advance measurements of individual characteristics that are related to biodynamic 
response can be used to predict surgical outcomes. In previous biomechanical simulations of 
15 refractive surgeries, the cornea was treated as a solid material. In this invention, the highly 
porous nature of and the normal negative pressure within the cornea is taken into 
consideration. By considering the changes in the position-dependent negative ("imbibition") 
pressure due to the tissue removal thorough surgery inside the cornea, the central portion of 
corneal surface after surgery is shown to flatten. Therefore, through the inventive method, 
20 which includes the position-dependent negative "imbibition" pressure, the basic clinically 
known results of LASIK or PRK can be correctly modeled. In addition, the clinically 
observed peripheral thickening after the LASIK or PRK is also demonstrated through the 
usage of simplified models. 

In the existing simulations for refractive surgeries, only the mechanical loading 
25 external to the cornea such as the intraocular pressure is considered. In the invented 
simulations, the cornea is basically modeled as a fluid-filled porous material as opposed to a 
solid material in the existing simulations. Therefore, the mechanical loading inside the 
cornea such as the negative ("imbibition") pressure applied on the internal pore boundary is 
directly incorporated into the simulations. Because of this major difference between the two 
30 methods, the inventive methodologies can reproduce and predict actual clinical physical 
phenomenon. 
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Other advantages and benefits of this invention will be apparent from the following 
detailed description. 

Drawings 

Figure 1 is a schematic of the simple "shape-substraction" paradigm for correction of 
5 myopia. 

Figure 2 is a schematic diagram of shape-subtraction model of refractive surgery for 
myopic ablation (left) and hyperopic ablation (right). 

Figure 3 presents post-operative minus pre-operative difference maps: elevation 
(upper left), pachymetric (lower left), tangential curvature (upper right) and axial (lower 
10 right). 

Figure 4 is a tangential curvature difference map (post-op minus pre-op) after a 
LASIK procedure. 

Figure 5 is a pachymetry difference map (post-op minus pre-op) after a myopic 
LASIK procedure. 

15 Figure 6 has post-op minus pre-op difference maps from a LASIK correction. 

Figure 7 has error maps of actual minus predicted post-op topography for the same 
subject as Figure 6. 

Figures 8 and 9 present maps of: A) pre-operative elevation (left) and tangential 
curvature (right); B) post-operative elevation (left) and tangential curvature (right); and C) 
20 elevation error (left) and tangential curvature error (right). 

Figure 10 presents regression analysis of peripheral stromal thickness of the superior 
region (upper plot) and inferior region (lower plot) against curvature. 

Figure 11 is a conceptual model that predicts biomechanical central flattening as a 
direct consequence of severed corneal lamellae. 

25 Figure 12 is a schematic diagram of the forces described by A) the interlamellar 

cohesive strength and B) the interlamellar shear strength. 

Figure 13 illustrates the impact of biomechanics on LASIK nomograms for primary 
and secondary hyperopia. 

Figure 14 shows a sample ablation depth profile calculated at 6 months after PRK. 
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Figure 15 is a comparison between planned ablation profiles and measured topographic 
elevation difference maps. 

Figure 16 has average difference maps between repeated measures of 20 eyes of 10 
normal subjects. 

5 Figure 17 has average difference maps between 1 day post-op LASIK and pre-op state. 

Figure 18 has average difference maps between 1 month post-op LASK. and pre-op 

state. 

Figure 19 presents composite difference maps of peripheral increases in curvature 
(upper right), elevation (upper left), and thickness (lower left) for patients who underwent 
1 0 LASIK surgery and had Orbscan corneal topography pre-operatively and at 6 months post- 
operatively. 

Figure 20 presents optical image views of the central region (A) and the inferior (B) 
superior (C) regions of the cornea. 

Figure 21 graphs central curvature shifts vs. regional peripheral stromal thickness 
15 changes. 

Figure 22 shows a pre-operative (upper left) and (lower left) curvature topography 
taken while the patient is supine for a LASIK procedure. 

Figure 23 the curvature difference map between pre-op and post-flap for the patient in 
Figure 22, calculated and displayed using customized software. The curvature difference map 
20 between post-flap and post-ablation for the same patient. 

Figure 24 is a model for an intact cornea, presented as a portion of an axisymmetric 
layered spherical shell fixed at the ends. 

Figure 25 is a model for cut cornea after PRK or LASIK surgery. 

Figure 26 illustrates flattening of the model cornea after surgery based on the fluid- 
25 filled porous material assumption of cornea. 

Figure 27 illustrates bulging of the model cornea after surgery based on the solid 
material assumption of cornea 

Figure 28 and 29 compare the overall corneal deformed shapes for intact and cut 
corneal models, with two different assumptions. 
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Detailed Description 

la order to improve the percentage of patients who achieve post-operative visual 
acuity of 20/20 or better, as well as allow minimization of aberrations, it is important to 
critically examine our current conceptual model of laser refractive surgery. There are three 
5 assumptions inherent in the shape-subtraction model that are not supported by the data and 
yet currently drive algorithm development. These flawed assumptions are: 1) the only 
portion of the cornea that is changed is within the ablation zone; 2) what you cut is what you 
get; and 3) even if there are changes outside the ablation zone, they don't affect central 
vision. 

10 Our comprehensive, integrated model of corneal response to laser refractive surgery 

takes into account three proposed components to the final corneal shape which determines 
vision: 1) ablation profile and laser parameters; 2) epithelial and stromal healing; and 3) 
biomechanical response to a change in structure. The third component, biomechanical 
response, has not been directly incorporated into existing ablation profiles in excimer laser 

15 refractive procedures. Therefore, this three pronged, systematic model, fully characterizing 
the corneal response represents an innovative approach to solving this complex problem. In 
addition, within the predictive model, an optimization approach is used to develop new 
ablation algorithms. Using optimization routines to produce the best possible outcome also 
represents a highly innovative approach to laser refractive surgery. This model and this 

20 optimization approach provide answers to the three assumptions or questions discussed 
above. 

Assumption #1 can be invalidated by anyone with access to a set of post-operative 
topography from PRK or LASIK. Outside the ablation zone, curvature significantly increases 
with the appearance of the characteristic red ring (high dioptric value) surrounding the central 

25 flattened zone after a myopic procedure. (Shown in the tangential map, not the axial.) In 
addition, elevation and pachymetry also increase outside the ablation zone, as measured by 
Orbscan topography, and shown with an example in Figure 3. Similar examples from 
another preliminary study are shown in Figure 4, which is a tangential curvature difference 
map (post-op minus pre-op) after a -12.5 diopter LASIK procedure using a Technolas 217 

30 (Munich, Germany) with a 5.5mm diameter ablation zone. The topography was acquired 
using an ORBSCAN I (Salt Lake City, Utah). The data were exported using the recorder 
function, and subsequently imported into custom software for analysis. Centrally, there is a 
decrease in curvature, as expected, indicated by the negative values and blue colors. The 

10 
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surrounding thin white area represents zero difference between the pre and post-op state. In 
the area outside the ablation zone, there is an unexpected increase in curvature which extends 
into the periphery, indicated by the positive values and red colors. Changes in corneal 
curvature clearly occur well beyond the ablation zone, challenging assumption #1. Figure 5 
5 is apachymetry difference map from the same patient. Inside the ablation zone, the thickness 
has decreased, as expected. However, outside the ablation zone, the thickness has 
unexpectedly increased. In addition, regression analysis between the central and peripheral 
curvature changes showed a significant negative correlation (p<0.0053), indicating that 
greater central flattening produced greater peripheral steepening. The peripheral increase in 

10 curvature is a known consequence of laser refractive surgery, but has never been fully 
explained other than as a "knee" at the edge of the ablation profile. If this were an accurate 
description, the change would be confined near the edge of the ablation zone, and yet it 
extends well beyond that region. Additional data in the form of elevation and pachymetry 
maps, from the same patient population, also showed significant peripheral increases in both 

15 elevation and pachymetry outside the ablation zone, corresponding to the increase in 
curvature. These paradoxical changes in elevation and pachymetry have been anecdotally 
attributed to measurement device errors, without exploration of the responsible 
biomechanical mechanisms. Figure 6 shows four difference maps from a sample patient in 
this population with a 6mm ablation zone, including elevation, tangential, axial, and 

20 pachymetry differences. The expected decreases in curvature, elevation and pachymetry 
centrally, within the ablation zone, are noted. However, the maps also show unexpected 
increases in elevation, pachymetry, and curvature outside the ablation zone. Therefore, the 
results of this study challenge the validity of the assumption that corneal curvature changes 
are confined to the ablation zone. The proposed mechanism for the paradoxical changes in 

25 elevation and pachymetry are the subject of this invention. 

The results of the analysis of assumption #1 demonstrated that shape changes occur well 
beyond the ablation zone. Assumption #2 can be examined more directly by comparing 
ablation profiles to actual corneal surface changes. In a preliminary study, 10 LAS1K 
patients without astigmatism who had symmetric ablation patterns, pre and post-operative 
30 topography was obtained using an Orbscan I. LASIK was performed using a Technolas 217 
excimer laser. Ablation algorithms were approximated using Munnerlyn's formulas. The 
calculated ablation profile was subtracted from the pre-operative topography to generate a 
predicted post-operative topography. Actual post-operative topography was compared to 

11 



SUBSTITUTE SHEET (RULE 26) 



WO 02/07660 




PCT/US01/22936 



predicted post-operative topography and error maps generated. RMS errors were calculated 
within the central 4mm diameter zone, outside the central zone, and over the entire map. 
Results averaged for all 10 subjects are given in Table 1, and demonstrate larger error outside 
the 4mm central zone than inside. 

5 TABLE 1 : RMS Error between Predicted and Measured Topography after LASIK 

elevation curvature 
mean (n=l 0) RMS error mean (n = 1 0) RMS error 

central 4mm diameter 18 ± 14 microns 4.92 ± 1 .89 diopters 

10 4-9mm diameter zone 23 ± 11 microns 8.06 + 1 .76 diopters 

overall 22 + 1 1 microns 6.85 + 1 .50 diopters 

The sample error map in Figure 7 shows a pattern of positive error peripherally in 
15 elevation, pachymetry, and curvature with negative error centrally in elevation, pachymetry 
and somewhat in curvature, although curvature is not consistently negative over the whole 
central region. These patterns of error cannot simply be attributed to using an estimated 
rather than known ablation profile, since they are both nonrandom and nonlinear. Ideally, the 
actual ablation algorithms should be used in a study of this type, rather than Munnerlyn's 
20 formulas, but they were not available due to their proprietary nature. However, despite using 
estimated ablation profiles, the evidence presented challenges the validity of assumption #2. 
These patterns of error are not consistent with a shape-subtraction model of ablation, but they 
are consistent with the biomechanical model to be proposed below. 

Assumption #2 can be further invalidated by examining the predicted topographic 
25 result with a known ablation profile, and comparing that to the measured result. In another 
preliminary study, Summit Technologies agreed to supply their proprietary ablation 
algorithms for several LASIK patients being treated with a Summit Apex Plus using a 
Krumeich-Barraquer microkeratome. The actual ablation algorithm was subtracted from the 
pre-operative topography (measured with an Orbscan II) to generate a ''predicted" post- 
30 operative topography. The predicted topography was subtracted from the measured post- 
operative topography to generate "error" maps, which are given in Figures 8 and 9 for two 
patients, along with the measured results. Several important features should be noted. First, 
the red area in the center of the tangential curvature error maps correspond to the central 
island treatment. This treatment involves delivery of extra pulses in the center of the cornea, 
35 which would cause excessive flattening IF they were not compensating for the central island 
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phenomenon - the genesis of which is not yet understood. Therefore, the curvature error 
maps appears to have a "central island" even though the post-operative topographies are 
smooth. Second, the red areas outside the ablation zone on the tangential curvature error 
maps correspond to the unexpected increases in curvature, which are not predicted by the 
5 shape-subtraction model, and are predicted via the biomechanical model presented herein. 

Finally, how can documented changes outside the ablation zone affect central 
curvature (assumption #3)? The answer to this question lies in the biomechanics of the 
corneal response to laser refractive surgery, which is unaccounted for in current ablation 
algorithms. It has been known since the inception of refractive surgery that altering the 

10 corneal structure will alter the shape of the entire cornea, whether using an incisional or 
thermal mechanism. Fundamentally, if the cornea were a piece of plastic, radial keratotomy 
would not have worked! Yet, with the development of laser refractive surgery, the structural 
link between the central and peripheral cornea was ignored. The cornea was thought of quite 
simply as a homogeneous structure to be "sculpted" into a new shape. However, this 

15 conceptualization cannot account for all of the corneal shape changes that occur after an 
ablative procedure. An important component of the proposed biomechanical model of 
corneal response to laser refractive surgery, to be described in the next section, is peripheral 
stromal thickening or an increase in corneal elevation outside the ablation zone. Evidence 
that this occurs was presented in the analysis of assumption #1. However, the question 

20 remains as to the linking of the central and peripheral corneal events. Is the peripheral 
increase in corneal elevation statistically correlated to central curvature changes, or are they 
independent, unrelated occurrences? Regression analysis between central curvature change 
and peripheral elevation change from the 30 subjects who underwent LASIK demonstrated a 
significant positive correlation (R 2 = 0.56, p < 0.0001), indicating that the greater the increase 

25 in elevation outside the ablation zone, the greater the curvature change (flattening) centrally. 
This provides evidence that the peripheral response of the cornea is linked to the central 
response, although it does not account for all the variance. 

The following case study further illustrates how central curvature changes can 
actually track peripheral thickness changes after ablation. The patient had 6mm diameter 
30 PRK (VISX Star) for a refractive error of -3.75 + 0.75 x 90, with Orbscan I topography pre 
and post-operatively. Both superior and inferior corneal thicknesses were calculated by 
averaging a 1mm diameter region in the pachymetry map along the vertical meridian, at the 
edges of a 7.5mm diameter circle centered on the ablation zone. Central curvature was 
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calculated by , averaging the curvature in the 3mm diameter central region of the tangential 
map. The values as a function of time relative to surgery are given in Table 2. 

TABLE 2: Case Study of Curvature and Thickness Changes after 6 mm Diameter 

PRK 

5 





Superior 


Inferior 


Central 


Central 




Thickness 


Thickness 


Thickness 


Curvature 




(microns) 


(microns) 


(microns) 




preop 


715 


682 


605.. 


41.55 diopters 


30 minutes 


770 


729 


710 


36.20 diopters 


day 4 


727 


691 


589 


39.39 diopters 


3 weeks 


744 


701 


568 


38.81 diopters 


3 months 


722 


690 


570 


39.32 diopters 



15 

Note that both central curvature and peripheral thickness demonstrate post-operative 
fluctuation in opposite directions. Regression analysis of central curvature vs peripheral 
thickness was performed, and the plots are given in Figure 10. Central curvature has a strong 
negative correlation with peripheral thickness, both inferior and superior, meaning the greater 
20 the peripheral thickness, the flatter the central curvature. This provides further support that 
the central and peripheral changes are linked and assumption #3 is not valid. 

Evidence has been presented challenging the inherent assumptions of the shape- 
subtraction model of refractive surgery. Clinical data indicate that substantial changes occur 
outside the ablation zone that are structurally linked to central curvature, and may affect 
25 central vision. The proposed mechanism for the measured increases in elevation, thickness, 
and curvature which occur outside the ablation zone will be presented below. 

A conceptual model is presented in FIGURE 11 that predicts biomechanical central 
flattening as a direct consequence of severed corneal lamellae. Rather than a piece of plastic, 
the cornea may be conceived as a series of stacked rubber bands (lamellae) with sponges 

30 between each layer (interlamellar spaces filled with ground substance or matrix). The rubber 
bands are in tension, since there is a force pushing on them from underneath (intraocular 
pressure), and the ends are held tightly by the limbus. The amount of water that each sponge 
can hold is determined by how tautly the rubber bands are pulled. The more they are pulled, 
the greater the tension each carries, the more water is squeezed out of the interleaving 

35 sponges, and the smaller the interlamellar spacing. This is analogous to the preoperative 
condition in Figure 11A. After myopic laser refractive surgery, a series of lamellae are 
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severed centrally and removed, as shown in Figure 11B. The remaining peripheral segments 
relax, just like the taut rubber bands would relax once cut. With the reduction of tension in 
the lamellae, the squeezing force on the matrix is reduced and the distance between lamellae 
expands (negative intra-stromal fluid pressure), analogous to the sponges taking up water if 
5 the rubber bands are cut. This allows the periphery of the cornea to thicken. Due to the 
crosslinking between lamellar layers, the expansion force pulls on the underlying intact 
lamellae, as indicated by the arrows pointing radially outward. An outward force in the 
periphery pulls laterally on the center and flattens it. Thus, the cornea will flatten centrally 
with any procedure that circumferentially severs lamellae. The biomechanical flattening 
10 enhances a myopic procedure, works against a hyperopic procedure, and will cause flattening 
in a "non-refractive' 9 PTK. This includes myopic profiles, hyperopic profiles, constant depth- 
PTK profiles, as well as the simple cutting of a LASIK flap. 

In every LASIK procedure, a flap is cut with a microkeratome to a thickness of 
approximately 160 microns. Biomechanically, this is very close to a 160 micron depth 
15 ablation, and according to the proposed model, should induce corneal flattening. The amount 
of flattening produced should be indicative of the biomechanical response of the cornea to 
refractive surgery. Topography of this epithelial surface of the cornea, before and just after 
cutting the flap, permits alteration of ablation algorithms in real time to account for individual 
biomechanics, and improves surgical outcomes. 

20 A Biomechanical Model of Keratectomy-Induced Curvature Change 

In recent decades, a great deal of effort has been devoted to characterizing the cornea 
as a biomechanical entity. This work has included in vitro measurement of key material 
properties such as the modulus of elasticity and shear modulus, and numerous computational 
models have been generated to approximate the structural response to simulated incisional 
25 keratectomy. While these efforts have undoubtedly advanced the basic understanding of 
corneal behavior under certain specific conditions, less has been done relative to modeling 
the structural changes induced by laser keratectomy, and the predictive value of existing 
numerical models is limited by their implicit simplifications (e.g. modeling the cornea as a 
solid substance). 

30 Mechanical Model of Keratectomy-Induced Corneal Flattening 

The basis for the proposed biomechanical theory of corneal response lies in the 
lamellar structure of the corneal stoma, which is altered in PTK, PRK, and LASIK. It is 
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known the stroma dominates the mechanical response of the cornea to injury. The stromal 
lamellae, by virtue of their extension across the corneal width and continuity with the corneo- 
scleral limbus, bear a tensile load arising from intraocular pressure and extraocular muscle 
tension. The propensity of the stroma to imbibe fluid and swell, which is attributed to the 
5 hydrophilic macromolecules of the extracellular matrix, is resisted by the corneal limiting 
layers, the metabolic activity of the corneal endothelium, and the compressive effects of 
lamellar tension. The thickness of the stroma is lineraly related to the hydration and thus is 
an accessible indicator of changes in corneal mechanical equilibrium. It is proposed that 
when tension-bearing lamellae are disrupted by central keratectomy, their peripheral 

10 segments relax, causing local decompression of the extracellular matrix, and a compensatory 
influx of stromal fluid adding to increases in interlamellar spacing and peripheral thickness. 
Increases in thickness are limited ultimately by interlamellar crosslinks, which are 
preferentially distributed in the anterior one-third and periphery of the stroma, as depicted 
schematically in Figure 11, and are postulated to contribute to regional differences in 

1 5 lamellar shearing strength and interlamellar cohesive strength. A pivotal aspect of the model 
is the prospect that interlamellar stress generated in the expanding stromal periphery is 
communicated through this lattice of crosslinks to lamellae comprising the postoperative 
optical surface. Arrows in Figure 11 indicate expansion of the peripheral stromal matrix 
secondary to ablation-induced lamellar relaxation. Peripheral thickening is proposed to exert 

20 an anterolateral tension at the ablation margin and stimulate central flattening. Coupling of 
peripheral and central tensile loads is mediated by interlamellar cohesive forces ("x") 
preferentially distributed in the anterior-peripheral stroma (shaded). Components of force in 
the direction of peripheral expansion result in central corneal flattening and peripheral 
steepening. 

25 During PRK, PTK and LASIK, central ablation causes an immediate circumferential 

severing of corneal lamellae under tension, with a subsequent relaxation of the corresponding 
peripheral lamellar segments. This causes a peripheral decompression of the extracellular 
matrix and an increase in stromal thickness outside of the ablation zone. It is further 
hypothesized that this response has important effects on central curvature due partly to the 

30 presence of interlamellar crosslinking, which predominates in the anterior one-third and 
periphery of the stroma and has been associated with significant lamellar shearing strength 
and interlamellar cohesive strength. Interlamellar stress is generated in the expanding 
stromal periphery and may be transferred through this network of crosslinks to the underlying 
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lamellae, whose central portions comprise the postoperative anterior surface (Figure 11). 
The outward expansive force in the periphery causes the central cornea to flatten, 
independent of the ablation profile. Thus, even in the absence of a myopic (centrally 
weighted) ablation pattern, biomechanical changes in the cornea may produce an acute 
5 postoperative flattening of the central cornea and a refractive shift toward hyperopia. This 
response is perhaps most clearly demonstrated by the clinical phenomenon of unintended 
hyperopic shift in PTK, but its contribution to refractive variability in PRK and LASHC is 
also important, since biomechanical central flattening will enhance a procedure to treat 
myopia and oppose a procedure to treat hyperopia. An additional feature of the model — an 
1 0 acute depth-dependent response — will also be discussed. 

Of the five anatomic layers of the cornea, (the epithelium, Bowman's layer, the 
stroma, Descemet's layer and the endothelium), only Bowman ? s zone and the stroma contain 
collagen fibrils. These layers are thus presumed to provide the majority of the corneas' 
tensile strength. Although tension in the superficial epithelial cells has been cited as a 
15 potential mechanism for maintaining a smooth optical surface in the presence of underlying 
stromal surface undulations, removal of the epithelium causes little or no change in the 
anterior corneal curvature, and the epithelium is generally attributed a minimal role in corneal 
tensile strength. 

The mechanics of Descemet's layer - the 7-um-thick hypertrophied basil lamina of 
20 the underlying endothelium — were studied in intact human and rabbit globes, and Compared 
to findings in stroma from the same species. Stress-strain investigations revealed that in the 
in situ human cornea, the stromal stress-strain curve is quite steep, corresponding to a high 
Young's modulus, while Descemet's layer is highly extensible. Superimposing the stress- 
strain curves for the two layers demonstrated that Descemet's layer is essentially unstrained 
25 over a large range of intraocular pressures for which the stroma is simultaneously fully 
strained. This leads to the conclusion that, like the epithelium, Decemet's layer probably does 
not bear a significant proportion of the corneal tension over a physiologic range of pressures. 
Marked differences in the stress- strain relationships of rabbit and human stroma highlight the 
need for critical interpretation when using the rabbit as a model for corneal mechanical 
30 behavior in humans. 

Structurally, Bowman's zone is little more than an acellular extension of the anterior 
stroma. Relative to the stroma, the individual collagen fibrils are two-thirds smaller (20 to 25 
nm in diameter) and more randomly oriented throughout the 8 to 12-um- thick sheet. 
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However, it has long been believed that Bowman's zone contributes a structural rigidity to the 
cornea that is distinct from that provided by the stroma. But recent studies, using a linear 
extensiometer to examine a number of mechanical parameters in deepithelialized corneal 
strips obtained from fresh human cadaver eyes have indicated that removal of Bowman's 
5 layer did not significantly alter the constitutive mechanical properties of the cornea. 



Finally, the stroma, which makes up about 90% of the total corneal thickness, is the 
layer thought to most significantly influence the mechanical response of the cornea to injury. 
The stroma is approximately 78% water by weight, 15% collagen, and 7% other proteins, 
proteoglycans, and salts. Three hundred to five hundred lamellae - flattened bundles of 

10 parallel collagen fibrils ~ run from limbus to limbus without interruption. In the posterior 
two-thirds of the stroma, the lamellae are successively stacked parallel to the corneal surface 
such that each lamellae has an angular offset from its anterior and posterior neighbors. 
Anteriorly, the lamellae are more randomly oriented, often obliquely to the corneal surface, 
are more branched, and are significantly interwoven. These regional differences correlate 

15 well with observations that shearing of the lamellae is generally difficult in the human 
stroma, but particularly so anteriorly. Collagen interweaving is also more extensive in the 
corneal periphery than in its center. 

It is important to distinguish between cohesive strength and interlamellar shear 
strength as the two terms are frequently used interchangeably. Cohesive strength has been 

20 measured as the force required to separate a stromal sample along a cleavage plane parallel to 
the lamellar axes by pulling in a direction perpendicular to the cleavage plane, like peeling a 
banana (as shown in Figure 12A). Thus, the cohesive strength is a measure of the 
interlamellar resistance to separation in the transverse direction and is expressed as a function 
of distance from the corneal center. Alternatively, the shear strength manifests as a 

25 resistance to shearing or sliding of one lamellae over another in the plane parallel to the 
lamellar axes (the longitudinal direction); as such, it is an integrated function of the 
connective forces across the entire lamellar interface and is therefore conceivably larger in 
magnitude (Figure 12B). Both forces are likely to contribute to the peripheral-to-central 
transfer of stress in the proposed biomechanical model of corneal response to laser ablation. 

30 Though the interlamellar cohesive strengths of the temporal and nasal peripheries of 

the horizontal meridian were equivalent in the study referred to above, cursory studies of the 
vertical meridian revealed large and consistent differences between the strengths of the 
superior and inferior regions. These differences were confirmed in a study which revealed a 
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mean central strength of 0.165 ± 0.0088 N/mm and mean peripheral strengths at 4-mm from 
the central cornea of 0.185 .± 0.0088 N/mm and 0.234 ± 0.0137 N/mm (p < .01) for the 
inferior and superior regions, respectively. Again, interlamellar cohesive strength in the 
human corneal stroma was shown to be greater in the periphery than in the center. 

5 These studies provided anatomical and mechanical evidence of an interlamellar 

cohesion strength at 50% depth that is greater peripherally than centrally and greater 
superiorly than inferiorly. The predominance of interlamellar binding in the anterior one- 
third of the stroma indicates that the magnitudes of these forces are even greater in that 
region. Whatever bearing these interlamellar relationships have on the shape of the 
10 structurally intact cornea, their importance is likely to increase considerably in the ablated 
cornea when a redistribution of tensile loading induces non-physiological stresses between 
cut and uncut lamellae as illustrated in Figure 11. Furthermore, the asymmetric distribution 
of cohesive strengths can be an important source of induced corneal astigmatism in 
keratorefractive procedures. 

15 The lamellar organization of the stroma and the capacity of this collagenous network 

to bear tension is a natural starting point for biomechanical models of curvature change. A 
relationship that is generally neglected in this context, however, is the relationship between 
the interfibrillary constituents of the stroma and water, the major component of the stroma. 
The collagen fibers are enmeshed in a ground matrix of glycosaminoglycans (GAG) such as 

20 keratan sulfate and chondroitin sulfates of varying degrees of sulfation. Both substances, but 
particularly chondroitin sulfate, are markedly hydrophilic and contribute to a negative 
intrastromal fluid pressure under which the entire stroma is heavily compressed. Intraocular 
pressure further compresses the stroma through its direct effect at the posterior surface and by 
its contribution to lamellar tension. The intrastromal pressure, often called the "swelling" 

25 pressure because of its tendency to draw water into the stromal ground substance, has been 
measured as -50 to -60 mmHg through a variety of in vitro and in vivo techniques. 

In the normal physiologic state, this swelling tendency is resisted and relative 
dehydration is maintained by a combination of lamellar tension, anterior evaporation of the 
tear film, low permeability of the epithelial and endothelial layers to water, and active 
30 endothelial transport of bicarbonate. During the act of central ablation, however, a number of 
lamellae proportional to the depth of ablation are obliterated centrally and tension is lost in 
their unablated peripheral segments. The resulting loss of compression introduces a 
hydrostatic disequilibrium, and the peripheral stroma thickens as it takes up fluid. While 
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tissue expansion is likely to occur within the ground substance rather than in the fibrils 
themselves, the source of the fluid has not been discerned. In situ, the cornea may imbibe 
additional water through the limbus as a result of the pressure difference between the 
perilimbai capillaries and the stromal interstitium. As one would expect, the swelling 
5 pressure diminishes as stromal hydration and thickness increase and a new steady-state is 
established. 

As stated earlier, the essential link between peripheral stromal thickening and central 
flattening is the existence of a mechanical relationship between disrupted and intact lamellae. 
Contrarily, if lamellae are assumed instead to be structurally and mechanically independent 

10 of neighboring lamellae — that is, arranged in layers without any interconnections — then the 
severed peripheral segments are unable to bear or transfer any tension. This is, in fact, a 
simplifying assumption incorporated into most numerical models of refractive surgery. 
Within this scenario, the tensile load previously borne by the full complement of lamellae is 
shifted to the remaining posterior fibers, which now strain (stretch) slightly under the 

15 concentrated stress. If the limbal circumference is assumed to be fixed, the stretch may not 
be accomplished by an increase in corneal diameter and must therefore occur as central 
corneal bulging and anterior steepening. However, if the proposed peripheral response is 
considered, this occurs as a peripheral thickening and peripheral steepening with coincident 
central flattening. Although the latter scenario aptly describes the characteristic peripheral 

20 "knee" and central flat zone of PTK, PRK and LASIK, finite element analyses of uniform 
thickness profiles have not incorporated the peripheral response and have predicted a corneal 
configuration opposite that produced clinically. In short, consideration of the proposed 
peripheral stromal response and its mechanical relationship to the central cornea, as well as 
the negative intrastromal pressure, is critical not only for correctly predicting the magnitude 

25 of refractive correction, but in some cases, the direction. 

Peripheral Thickening Affects Central Curvature 

One of the important implications of this work is the notion that shape changes 
measured outside of the ablation zone can have a significant impact on curvature changes 
within the ablation zone. The ablation zone analysis demonstrated in a controlled fashion that 
30 central curvature change in refractive surgery is not solely a product of the ablation pattern. 
When peripheral thickness and ablation zone bias were both included in a regression model, 
over 83% of the variance in curvature response accounted for by the model was explained by 
peripheral thickening. Although the impact of ablation pattern on acute curvature change is 
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likely to become increasingly influential when larger dioptric corrections are attempted, the 
associated increases in ablation depth would also be expected to further exacerbate the 
biomechanical response. Consequently, the utility of a pure shape-subtraction model is 
probably limited to cases involving only superficial ablation (i.e., correction of low refractive 
5 error or astigmatism in PRK, and removal of subepithelial scars in PTK). Finally, the argument 
for the role of lamellar interconnections playing an important role was bolstered by the 
correlations between central flattening and thickness changes in individual peripheral 
subregions. These correlations were strongest in the superior periphery, weakest in the inferior 
periphery, and strongest overall in the far superior periphery — a pattern that consistently 
1 0 reflects the strength distribution of interlamellar cohesive forces reviewed above. 

Hyperopic vs Myopic Corrections 

The proposed biomechanical model predicts additional flattening over and above 
whatever ablation profile is programmed, whether myopic with an intent to flatten, hyperopic 
with an attempt to steepen, or nonrefractive PTK. This would theoretically make hyperopia a 

15 more difficult procedure, since the biomechanical flattening would be in opposition to the 
ablation profile, thus requiring a deeper ablation to offset the biomechanical response,. This 
prediction was verified in the first 8 patients treated with the Autonomous Custom Cornea, 
wavefront-guided ablation procedure. Of the 8 patients treated, 5 were myopic and 3 were 
hyperopic. Only the wavefront data were used to program the laser. Neither refraction, nor 

20 empirical experience were used. All 5 of the myopic patients were overcorrected (additional 
biomechanical flattening in combination with programmed ablative flattening) and all 3 
hyperopic patients were undercorrected (biomechanical flattening in opposition to 
programmed ablative steepening). 

In a comparison of primary vs secondary hyperopic treatments (overcorrected post- 
25 myopic procedures), primary hyperopic procedures require a treatment of up to 35% greater 
. than the spherical equivalent, meaning greater than expected depth must be ablated to achieve 
the desired correction. This is consistent with the proposed model predicting biomechanical 
flattening, independent of the ablation profile. Greater depth must be ablated with a primary 
hyperopic treatment to overcome the biomechanical flattening which creates a surface shape 
3 0 effect opposing that of the ablation profile. 

The secondary hyperopic group required substantially less depth of ablation to 
achieve the same level of correction as the primary hyperopia group. This is also completely 
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consistent with the proposed biomechanical model. The secondary post-LAS IK hyperopic 
group already had an altered corneal structure with an associated biomechanical response 
from the first refractive procedure. In other words, there was less biomechanical flattening to 
overcome in the second procedure, thus requiring less depth of ablation for the same 
5 correction, since the biomechanical effect had already occurred in conjunction with the first 
procedure, as illustrated in Figure 13. Third, the secondary hyperopes had a more stable 
longer term refraction, with the primary hyperopes exhibiting slight regression over the 6 
month post-operative follow-up period- One possible explanation for this difference lies in 
the long-term healing related to the biomechanical effect, which would be exaggerated in the 
10 primary hyperopia group, as opposed to the secondary hyperopia group where the 
biomechanical effects had already occurred and were thus reduced. 

Peripheral Thickening After Laser Refractive Surgery In Vivo 
Indirect evidence of peripheral thickening after laser refractive surgery has been 
reported, which is consistent with the proposed model of biomechanical response. Munger, 

15 et. al. presented clinical results using a VISX Star system on a total of 25 eyes from 25 
patients with low astigmatism, treated with myopic PRK. (Munger R, Jackson W B, 
Mintsioulis G, Ablation Profile and Epithelial Regrowth after Myopic PRK with the VISX 
Star. American Society of Cataract and Refractive Surgery Annual Meeting, 1999) Corneal 
maps were acquired with PAR Corneal Topography System (CTS) at 2 hours before surgery; 

20 2-3 minutes after surgery (Sx); on the day the bandage contact lens was removed, and 1 
week, 3 months, and 6 months post-operatively. Ablation depth was defined as preop - post- 
op, and epithelial depth was defined as Sx - post-op, recognizing that edema may have been a 
factor in the immediate post-op map. Munger reported epithelial thickening at the ablation 
edge. However, an increase in the calculated epithelial measurement outside the ablation 

25 zone may be due to stromal thickening rather than epithelial, since the PAR CTS measures 
surface elevation and cannot distinguish between epithelial and stromal thickening. In 
addition, more convincing evidence lay in the ablation depth measurements as a function of 
distance from the apex at 6 months, an example of which is shown in Figure 14. Ablation 
depth is a more robust measurement than the much smaller epithelial thickness, and by 6 

30 months the cornea should be relatively stable. Note that outside the ablation zone, there is 
negative ablation depth, or increased elevation compared to the pre-operative state. This 
provides indirect evidence of peripheral thickening, an important feature of the proposed 
biomechanical model, with a different laser/topographer combination. 
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Sborgia, et al, has reported results of a corneal topography guided system called 
Corneal Interactive Programmed Topographic Ablation (CEPTA). (Sborgia C, Alessio G, 
Boscia F, Vetrugno M. Corneal Interactived Programmed Topographic Ablation: 
Preliminary Results. American Society of Cataract and Refractive Surgery Annual Meeting 
5 1999) A total of 28 eyes of 28 patients with irregular astigmatism were treated with a 
Laserscan 2000 and had pre and post-operative topography with an Orbscan. Two examples 
of the planned ablation profile and actual elevation difference maps are given in Figure 15. 
In both subjects, the left image is the planned ablation pattern, and the right image is the 
actual elevation difference map. Note for both cases the large green area outside the ablation 
10 zone in the maps on the left, where the cornea will not be ablated. However, the elevation 
difference maps on the right demonstrate not only a decrease in central elevation (shown in 
blue) which matches the planned profile, but are also accompanied by an unexpected increase 
in elevation outside the ablation zone (shown in red and yellow), again as predicted by the 
proposed biomechanical model. 

1 5 Characterization of the Biomechanical Response 

In a study to characterize the biomechanical corneal response to LASIK, and separate 
the resulting refractive effect from that produced by the ablation profile, 8 eyes of 4 subjects 
have been enrolled and analyzed. LASIK is performed with a Summit Apex Plus using a 
Krumeich-Barraquer microkeratome to make the flap. Corneal measurements are acquired 

20 pre and postoperatively with optical coherence tomography and the following corneal 
topographers for validation of effect: EyeSys, Humphrey Atlas, Keratron, Orbscan II, PAR, 
Technomed C-Scan, and TMS-L Average pre-operative refractive error of the small 
preliminary group is -6.875±2.03 diopters sphere + 0.8125+0.51 diopters cylinder. Only the 
Orbscan corneal topography data will be presented here. A comparison control group of 20 

25 eyes of 10 subjects had repeated Orbscan I topography acquired at intervals of from 1-2 days 
to establish baseline variations. The anterior tangential, anterior elevation and pachymetry 
data were exported to a customized topography tool for analysis and the cornea was divided 
into three regions: central 2.75 mm radius (5.5mm diameter), transition zone from a radius of 
2.75 - 3.25mm (5.5-6.5 mm diameter), and outside the ablation zone from a radius of 3.25 - 

30 4.5mm (6.5-9.0mm diameter). The pre-operative topography was subtracted from the post- 
operative topography for the surgical patients, and the repeated measurements were 
subtracted for the normal subjects. For the elevation maps, the two surfaces were fit within 
the 0.5mm transition zone. For all maps, average regional differences were calculated over 
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the normal and surgical populations, and statistical analysis was performed using the 
ANOVA ("Analysis of Variance") procedure in the software package, (Statistical Analysis 
System, Cary, NC) . 

Figures 16 and 17 show the composite difference maps of all subjects for the normal 
5 and 1 day post-operative surgical groups, respectively. Significant (p< 0.05) decreases in 
elevation, pachymetry, and curvature in the central zone between normals and surgical 
subjects were demonstrated. In addition, significant (p< 0.05) increases in elevation, 
pachymetry, and curvature in the outer zone were found, as predicted by the proposed 
biomechanical model. These increases persisted at 1 month post-operatively, as seen in 
10 Figure 18. 

A large retrospective study was conducted of 2,380 patients who received LASIK 
using a Technolas 217 laser, who had Orbscan Corneal Topography pre-operatively and at 6 
months post-operatively. Figure 19 shows the composite difference maps generated from 
the Orbscan data, demonstrating the same pattern of peripheral increases in curvature, 
15 elevation, and thickness over a large population. In addition, central curvature difference was 
highly correlated with peripheral elevation difference in a regression analysis (R 2 = 0.76, p< 
0.0001), consistent with the proposed biomechanical model of corneal response. 

In Vitro Studies 

Two in vitro studies were conducted. In the first study, 14 deepithelialized eye bank 
20 globes from 7 donors were subjected in paired-control fashion to either broad-beam PTK 
(-100-um depth, no programmed dioptric change) or sham photoablation. Changes in 
anterior curvature were measured by autokeratometry. Changes in stromal thickness in the 
vertical meridian were measured using corneal optical section image analysis. This is a 
technique using a modified endothelial camera to obtain cross-sectional corneal images, an 
25 example of which is shown in Figure 20. The corneal cross-section was divided into 5 
regions for analysis, far superior, near superior, central, near inferior, and far inferior. 
Analysis included evaluating peripheral thickness changes and geometric bias as predictors of 
curvature change. Geometric bias was defined as either a myopic or hyperopic bias in the 
pattern of ablation zone thickness loss in order to investigate current shape-subtraction 
30 theories of hyperopic shift in PTK. 

Results of this first study demonstrated that photoablation caused significant 
reductions in keratometric curvature (-6.28 ±3.23 D, P = .002) relative to untreated paired 
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controls. The mean keratometric shift measured during sham PTK was not significantly 
different from zero (+0.31 ± 0.85 D, P = .38). In addition to dramatic flattening of the 
spherical component of curvature, ablated corneas demonstrated significantly higher absolute 
magnitudes of keratometric cylinder (2.98 ± 0.88 D) than controls (0.46 ± 1.72 D, P - .009), 
5 indicating a potential biomechanical component to induced cylinder, potentially due to 
nonuniform distribution of corneal crosslinking. The relative peripheral stromal thickness 
change, expressed as the mean pairwise difference (PTK-control), was +57.3 ±42.8 um (P = 
.01) or +8.5 ± 5.7 % (P = .01), demonstrating significant thickening relative to controls and 
supporting the proposed new theory. Central curvature shifts were linearly dependent upon 

10 regional peripheral stromal thickness changes in ablated and control eyes (Figure 21). 
Correlations were significant in 3 of the 4 peripheral subregions and marginally significant in 
the far inferior subregion (P = .05). The far superior stromal subregion exhibited the 
strongest correlation to hyperopic shift (r = -0.70). Correlation coefficients are negative and 
indicate a tendency for. corneas to flatten more extensively with increasingly positive 

15 peripheral thickness changes, a result that is consistent with the proposed role of peripheral 
thickening as a mechanical stimulus for central flattening. The correlation between 
keratometric shift and the mean total peripheral stromal thickness change (the mean change 
over all 4 subregions) was significant (r = -0.67, P - .01), and equivalent results were 
observed when thickness changes were expressed as percent changes. 

20 On the other hand, calculated geometric bias favored flattening in only 50% of cases, 

despite the induction of hyperopic shift in all ablated corneas. In addition, mean geometric 
biases in ablated corneas and paired controls were not significantly different from zero (P = 
0.66 and 0.52, respectively) or from each other (P=0.74). In sharp contrast to the apparent 
relationship between increasing peripheral stromal thickness and decreasing central 

25 curvature, geometric bias demonstrated an absence of any linear relationship with anterior 
corneal flattening (r = 0.06, P = .85), Multiple transformations of the data were investigated 
and yielded no demonstrable correlations, challenging current shape-subtraction theories of 
hyperopic shift in PTK. The bottom line was that geometric bias was a poor lone predictor of 
anterior corneal flattening. The results of this study demonstrate unequivocally that epithelial 

30 hyperplasia and stromal remodeling are not absolute requirements for the production of 
clinical magnitudes of anterior flattening in PTK. 

A second in vitro study was performed to further examine the relationship between 
peripheral stromal thickening and central flattening. Specifically, a paired-control human 
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donor eye study (n=20) was conducted to assess the efficacy of preoperative topical 
glutaraldehyde (GTA) treatment as a technique for inhibiting PTK-induced peripheral stromal 
thickening and, secondarily, for attenuating the acute corneal flattening response. Each eye 
was individually mounted in a custom holder, inflated to normal intraocular pressure (15 
5 mmHg) and deepithelialized. According to a crosslinking protocol developed in preliminary 
experiments, one cornea of a given donor pair was immersed in a 15% dextran solution for 40 
minutes then transferred to 4% GTA/dextran for an additional 20 minutes; the fellow control 
was exposed to 15% dextran for 60 minutes. Each eye was subsequently subjected to 1) 
sham PTK, a same-eye control phase incorporated to account for thinning due to 

10 intraoperative dehydration, 2) PTK (5-mm-diameter, lOOum-depth) and 3) a 1-hour 
hypo-osmotic soak phase designed to assess the anti-swelling activity of stromal crosslinking. 
A scanning-slit topography system (Orbscan) was used to acquire triplicate thickness and 
curvature measurements before and after each experimental phase. Crosslinking significantly 
inhibited peripheral stromal thickening during PTK and postoperative hypo-osmotic 

15 immersion. In addition, during PTK, crosslinked corneas demonstrated 36% less hyperopic 
shift relative to paired controls (p = .001). The magnitude of this latter effect was linearly 
dependent upon the magnitude of crosslink-mediated suppression of the peripheral thickening 
response to PTK (r = 0.68, p =.03). The results demonstrate that acute hyperopic shifts in a 
donor model of PTK can be significantly reduced through preoperative application of a 

20 collagen crosslinking reagent and therefore support the conclusion that mechanical events in 
the corneal periphery play an important role in keratectomy — induced central curvature 
changes. This mechanism of curvature change is not accounted for in current surgical 
algorithms, and potentially plays a role in the variability in outcomes. 

Optical Coherence Tomography: 

25 Optical coherence tomography (OCT), is an imaging technology which has only 

appeared in the literature since the 1990s. OCT images are analogous to ultrasound images, 
in that they are two dimensional cross sectional images depicting tissue reflectivity. 
However, in the case of OCT, the images show infrared reflectivity and demonstrate a 10 
micron longitudinal resolution, as compared to the 50 micron longitudinal resolution 

30 associated with traditional ultrasound. While high frequency ultrasound does allow a 20 
micron longitudinal resolution, it cannot penetrate more than 4 mm beyond the cornea. The 
OCT also is also limited to 3 mm penetration. However, the optically clear nature of the eye 
allows the retina to be imaged to a depth of 3mm. For corneal imaging, the OCT offers 
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additional advantages in that it is non-contact, unlike high frequency ultrasound which 
requires a water bath and topical anesthetic. 

The principles by which OCT works are those of interferometry. The OCT is 
essentially a Michelson interferometer with the subject's eye placed at the end of one of the 
5 light paths. 

Prediction of Biomechanical Response via 
Corneal Pertubation and Biomechanical Modeling 

Most recently, in cooperation with two surgeons, we have measured the shape of the 
LASIK flap using a Keratron Scout, which is a portable, Placido-based topographer. It can 

10 be used in a vertical position intra-operatively to acquire the topography of the flap 
immediately after it is cut, before it is reflected, and thus prior to ablation. Figure 22 
demonstrates that the characteristic "red ring" on this patient is a biomechanical 
phenomenon, since it appears BEFORE any tissue is removed. Cutting the flap alters the 
corneal structure in the same way described by the model presented. The difference is that 

15 the severed lamellae are not ablated, but put back in place. Therefore, the stability of the 
greater than 1 diopter average central decrease in curvature shown by the difference map in 
Figure 23, which is generated simply by cutting the flap, is not known. Investigations are 
currently underway. 

Corneal measurements following the cutting of the corneal flap (or other 
20 perturbation), but before ablation, are a key factor in embodiements of this invention. As 
demonstrated above, the microkeratomic incision for the flap produces definite changes in the 
cornea, regardless of any subsequent removal of tissue. The redistribution of strain caused by 
the keratomic incision causes the central cornea to flatten and the peripheral stromal matrix to 
thicken and become steeper. This reshaping assists with a myopic correction, i.e., a 
25 correction where increased corneal curvature is prescribed, and works against a hyperopic 
correction. 

In one embodidment of practicing this invention, corneal topography, optical 
coherence tomography, ultrasound, refraction and/or wave front analysis is used both before 
and after the microkeratomic incision for the corneal flap (or other corneal perturbation), 
30 from a statistically sufficient number of subjects who may be surgical patients, and the 
differences compared to expected and achieved post-operative results (undercorrected or 
overcorrected). Once developed, comparison of the measurements before and after the 
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incision will allow one to adjust the ablative procedure to account for the predicted 
biomechanical response, as a function of the changes measured after cutting the flap (or other 
corenal perturbation).. Preferably, these adjustments are made by laser manufacturers to their 
laser algorithms. However, since existing laser algorithms do not take the biomechanical 
5 changes into account, this invention may also be practiced by developers of corrective tables 
for current laser systems, or by surgeons using such tables. Ablation profile adjustments can 
be made in advance of the ablation in a separate procedure, or in real-time as an intra- 
operative adjustment, after the pertubation, but before the ablation. 

As mentioned above, optimal surgical procedures also require consideration of the 
10 biomechanical results of the ablative procedure itself. This is best performed by: 

(1) taking preoperative measurements; 

(2) talcing post-pertubation measurements after the microkeratome cut (or other 
, corneal perturbation), but before the ablative procedure; 

(3) comparing the difference between the pre-operative and post-perturbation 
15 measurements; and 

(4) making appropriate adjustments in the ablative algorithm based on the 
magnitude and pattern of the post-pertubation response; and 

(5) taking additional measurements after the ablative procedure to document the 
results. La a second embodiment of practicing this invention, corneal topography, optical 

20 coherence tomography, ultrasound, refraction, and/or wave front analysis are acquired both 
pre-operatively and post-operatively, after the laser refractive procedure, from a statistically 
sufficient number of subjects who are surgical patients. The pre-operative data are used to 
partially define the specifications of individual biomechanical mathematical models of 
individual corneas (thickness profile, curvature profile, corneal size). The ablation profile 
25 with which an individual was treated, is then used to mathematically "remove" layers from 
the model, and the post-operative data are used to define the final condition of the model. 
PRK, LASIK, and LASEK would each require distinct mathematical models, since the tissue 
removed is distinctly located within the corneal stroma, (e.g. PRK and LASEK are surface 
ablations and LAS1K is a deeper ablation.) Knowledge of the pre-operative and post- 
30 operative state, measured from the subject, allows the model properties (e.g. Young's 
modulus) to be adjusted in order to match the final predicted post-operative state achieved by 
the model with the actual measured post-operative state, in an iterative process. The material 
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properties thus determined for all corneas and all models within the subject population will be 
correlated with the pre-operative population data (e.g. age, sex, race, years of contact lens 
wear) and measured pre-operative data (e.g. thickness, curvature, wavefront, corneal size), to 
determine which pre-operative parameters (both characteristic and measured) best predict the 
5 material properties and thus response. These correlations will be used to compile a program 
that will produce material properties as an output, with pre-operative data as input. Once 
developed, the material properties program, along with individual pre-operative data, will be 
used to fully define a predictive biomechanical model of individual response to corneal 
ablation. Model predictions will allow one to adjust the ablative procedure to account for the 
10 predicted biomechanical response. Preferably, these adjustments are made by laser 
manufacturers to their laser algorithms. However, since existing laser algorithms do not take 
the biomechanical changes into account, this invention may also be practiced by developers 
of corrective tables for current laser systems, or by surgeons using such tables. 

As mentioned above, optimal surgical procedures also require consideration of the 
1 5 biomechanical results of the ablative procedure itself. This is best performed by: 

(1) taking preoperative measurements; 

(2) predicting the biomechanical properties based on pre-operative data; 

(3) fully defining a predictive biomechanical model based on pre-operative data 
and predicted material properties; and 

20 (4) making appropriate adjustments in the ablative algorithm based on model 

predictions; and 

(5) taking additional measurements after the ablative procedure to document the 
results. 

In a third embodiment of practicing this invention, corneal topography, optical 

25 coherence tomography, ultrasound, refraction, and/or wave front analysis are acquired pre- 

operatively, both before and after the microkeratomic incision for the corneal flap (or other 

corneal perturbation), and post-operatively, after the laser refractive procedure, from a 

statistically sufficient number of subjects who are surgical patients. The pre-operative data 

and are used to partially define the specifications of individual biomechanical mathematical 

30 models of individual corneas (thickness profile, curvature profile, corneal size). The ablation 

profile with which an individual was treated, is then used to mathematically "remove" layers 

from the model, and the post-operative data are used to define the final condition of the 

29 



SUBSTITUTE SHEET (RULE 26) 



WO 02/07660 




PCT7US01/22936 



model PRK, LASIK, and LASEK would each require distinct mathematical models, since 
the tissue removed is distinctly located within the corneal stroma, (e.g. PRK and LASEK are 
surface ablations and LASIK is a deeper ablation.) Knowledge of the pre-operative, and 
post-operative state, measured from the subject, allows the model properties (e.g. Young's 
5 modulus) to be adjusted in order to match the final predicted post-operative state achieved by 
the model with the actual measured post-operative state, in an iterative process. The material 
properties thus determined for all corneas and all models within the subject population will be 
correlated with the pre-operative population data (e.g. age, sex, race, years of contact lens 
wear), measured pre-operative data (e.g. thickness, curvature, wavefront, corneal size), as 

10 well as the post-pertubation data, to determine which pre-operative parameters (both 
characteristic and measured) and/or post-pertubation parameters best predict the material 
properties and thus response. These correlations will be used to compile a program that will 
produce material properties as an output, with pre-operative and post-pertubation data as 
input. Once developed, the material properties program, along with individual pre-operative 

1 5 data and post-pertubation data, will be used to fully define a predictive biomechanical model 
of individual response to corneal ablation. Model predictions will allow one to adjust the 
ablative procedure to account for the predicted biomechanical response, based on pre- 
operative and post-pertubation data. Preferably, these adjustments are made by laser 
manufacturers to their laser algorithms. However, since existing laser algorithms do not take 

20 the biomechanical changes into account, this invention may also be practiced by developers 
of corrective tables for current laser systems, or by surgeons using such tables. Ablation 
profile adjustments can be made in advance of the ablation in a separate procedure, or in real- 
time as an intra-operative adjustment, after the pertubation, but before the ablation. 

As mentioned above, optimal surgical procedures also require consideration of the 
25 biomechanical results of the ablative procedure itself. This is best performed by: 

(1) taking preoperative measurements; 

(2) taking post-pertubation measurements after the microkeratome cut (or other 
corneal perturbation), but before the ablative procedure; 

(3) predicting the biomechanical properties based on pre-operative data and post- 
3 0 pertubation data; 

(3) fully defining a predictive biomechanical model based on predicted material 
properties, pre-operative data; and post-pertubation data 
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(4) making appropriate adjustments in the ablative algorithm based on model 
predictions; and 

(5) taking additional measurements after the ablative procedure to document the 
results. 

5 These processes may be performed most effectively with laser manufacturers. They 

have already made empirical adjustments to the algorithms originally proposed by Munnerlyn 
et al, and should be eager to make further adjustments that would produce superior results. 
Similarly, the results of the corneal healing process should be considered. This is another 
iterative process taking measurements immediately following an operation and at specified 
10 intervals thereafter to determine biomechanical modifications in the cornea following 
surgery. This is also applicable to customized ablation, based on topography-guided or 
" wavefiront-guided procedures that do not rely on the Munnerlyn formulas. 

The corneal healing response is characterized using OCT imaging to measure 
epithelial, stromal, and flap thickness. For each examination, five images are acquired in the 

15 vertical meridian and five in the horizontal meridian, until the remote site completes the 
software to scan the cornea in a radial spoke pattern covering the central cornea within an 
8mm diameter. The OCT examinations is performed at the same time points as the 
topographic measurements for comparison purposes: pre-operatively and post-operatively at 
one day, one week, one month, three months, and six months, with the addition of one time 

20 point immediately after surgery in order to visualize the flap interface. 

The biomechanical response is believed to be immediate, with modification through 
the healing phase. Therefore, the behavior of the model is defined based on one day, one 
week, one month, 3 months, and 6 months post-operative measurements, with the exception 
of the OCT measurements that will be taken immediately after the surgery in order to 

25 visualize the flap. Once the model is fully characterized, and performance is satisfactory, an 
optimization approach is used to design new ablation algorithms. Rather than providing an 
ablation profile to the model and predicting the corneal response, optimization criteria is used 
as input to produce the necessary ablation profile to meet the criteria. Examples may be to 
maximize ablation zone diameter, while minimizing ablation depth to reach a target central 

30 curvature with minimal aberrations. 

Post-operative corneal shape, and thus visual performance, is the function of at least 
three factors: the ablation profile, the healing process, and the biomechanical response of the 
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cornea to a change in structure. Only by increasing our knowledge of the interaction of these 
factors can predictability in PRK and LASIK be improved. This has important implications 
in the development of new ablation algorithms and guided procedures. It points to an 
optimization approach, rather than a priori defining an "ideal" corneal shape that is ultimately 
5 not achievable. There are only certain shapes a cornea will biomechanically assume. For 
example, the deeper the peripheral cut in a myopic procedure to generate a potentially 
desirable post-operative "prolate" shape, the greater the number of severed lamellae and the 
greater the biomechanical central flattening response to counter the effect. Both the ablation 
profile and the biomechanical response need to be taken into account, as well as the healing 

10 response. Therefore, step one is to gain a better understanding of corneal response to 
standard ablation profiles in well-controlled studies, before moving into the realm of 
customized procedures. The outcome measures in these controlled studies should be more 
comprehensive than in the past to allow us to thoroughly interrogate the corneal response. 
This means we need to measure and report the outcome eiror in terms of the predicted 

1 5 topography and/or the predicted wavefront, not just visual acuity, spheres and cylinders. 

Sample Finite Element models demonstrating the unique attributes of the invention 

Major advantages of this invention are obtained by directly introducing the negative 
("imbibition") pressure into the simulations with simple models as shown in Figures 24 and 
25. In Figure 24, a model for intact cornea is presented as a portion of an axisymmetric 
20 layered spherical shell fixed at the ends. This model consists of alternately placed four thin 
hard layers (layer #1,3,5,7) and three thick soft layers (layer #2,4,6). In this simplified 
model, the collagen fibers are schematically modeled as impermeable hard layers. Soft 
layers, which are assumed to be highly porous and fully saturated with water, representing the 
matrix or ground substance containing proteoglycans. 

25 As discussed before, the major portion of the actual cornea (stroma) consists of over 

300 ~ 500 layers throughout the thickness and each layer is reinforced by a number of long 
collagen fibers with the same direction. The direction of the collagen fibers in each layer is 
different so that the in-plane reinforcement is more or less uniform. The model is simplified 
not only because the number of layers is reduced from 300-500 to 7 but also because the 

30 collagen fibers are collectively assumed to be thin hard shell layers instead of maintaining 
individual long cylindrical shape. This simplified model, however, is significantly different 
from any existing corneal finite element model in the sense that the basic microstructure of 
the cornea as a layered or lamellar structure, along with being highly porous, is directly 
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incorporated into the finite element model. In addition, basic substances which make up the 
cornea, such as collagen fibers and ground substance (proteoglycans) and water are modeled 
separately. In Figure 25, a simple model for an ablated cornea after surgery relevant to PRK 
is presented. A part of the two outer layers (layer #1 and #2) is removed from the intact 
5 corneal model. 



TABLE 3: Layer Geometry and Material Properties 



Layer number 


#1,3,5,7 
(Hard) 


#2,4,6 
(Soft) 


Thickness [mm] 


0.02 each 
(4 layers) 


0.15 each 
(3 layers) 


Young's modulus * 
EHPal 


5xl0 8 


4xl0 4 


Poisson's ratio ** 

V 


0.3 


0.3 



* J.O. Hjortdal, Regional Elastic Performance of the Human Cornea, Journal of 
Biomechanics (1996) 29, 931-942. 

10 **TJ. Shin, R.P. Vito, L.W. Johnson and B.E. McCarey, The Distribution of Strain in the 
Human Cornea (1997) 30, 497-503. 

In the above simple models, the overall material properties for both hard and soft 
layers are assumed to be linearly elastic and isotropic. The elastic moduli are taken from the 

15 literature. As discussed before, the soft layer is assumed to be a highly porous material. 
Based on the existing theories for fluid-filled porous materials (Katsube, 1985; Katsube and 
Carroll, 1987a,b), this porous material assumption requires additional elastic material 
parameter such as the bulk modulus of the solid matrix material of the soft layer without 
pores. The physical measurement of the bulk modulus of the matrix itself (without pores) 

20 will be extremely difficult and has not been reported in literature. In the simple models, this 
bulk modulus is assumed to be infinite, and the matrix itself (without pores) is assumed to be 
incompressible. This assumption is reasonable because the volume fraction of the matrix is 
extremely low in the cornea, and much of the volume compaction or expansion of the soft 
layers can be attributed to the pore space compaction or expansion through the deformation 

25 of the matrix. 

The volumetric expansion of a porous solid material due to the internal pore pressure 

is similar to volumetric thermal expansion due to temperature increase although the physical 
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mechanisms are completely different. Therefore, this volume change due to the internal pore 
pressure can be replaced by the thermal expansion term in any of the commercially available 
finite element codes, provided the material properties of the (dry) porous material (without 
fluid) and those of the solid matrix (without pores) are known. 

5 As shown in Figure 24, the negative ("imbibition") pressure of magnitude 60 mm Hg 

is applied uniformly throughout the pore boundary of highly porous layers #2, #4, and #6, in 
addition to the intraocular pressure of magnitude 15 mmHg applied on the posterior surface 
of the cornea. In carrying out the simulations, the commercially available finite element code 
ABAQUS with four-node axisymmetric elements is employed. In Figure NK2, the negative 

10 ("imbibition") pressure is applied only to porous layer #4 and #6, and zero pore pressure is 
assumed in the centrally ablated or cut porous layer #2. This simplified assumption is based 
on the fact that the fluid movement is most likely to occur in the in-plane direction rather than 
out-of-plane direction due to the layered nature of actual cornea. Since the internal fluid is 
suddenly exposed to the atmospheric pressure at the place where the layer is severed, the 

1 5 internal pressure is set equal to zero in layer #2. 

Results based on the sample finite element models 

The contours of the spherical surface B (the interface between the soft layer #2 and 
the hard layer #3) for the intact cornea model (Figure 24) and ablated or cut cornea model 
(Figure 25) are plotted in Figure 26. The displacements are magnified by a factor of 10 so 
20 that the comparison before and after the surgery can be visually observed. The spherical 
surface B flattens after the surgery. This flattening due to the surgery demonstrates the same 
trend as the conceptual model, supported by clinically observed results described earlier in 
this application. 

As discussed previously, in all the existing methodologies of refractive surgeries, the 
25 fact that the cornea is a porous material is ignored. Instead, the cornea is treated as a solid 
material. Therefore, the negative ("imbibition") pressure applied through the internal pore 
boundary inside the cornea is not taken into account. In order to elucidate these differences, 
we plot the deformation of the spherical surface B before and after the surgery by setting 
negative "imbibition" pressure in soft layers #2, #4, and #6 equal to zero. As shown clearly 
30 in Figure 27, the spherical surface B bulges after the surgery. This is not consistent with 
clinical results of surgery. The bulging occurs because the central part of the cornea has less 
rigidity due to the removal of tissue and therefore deforms more after the surgery. 
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In order to demonstrate further the differences between existing simulations and 
simulations based on the invention, the overall corneal deformed shapes before and after the 
surgery are compared based on the two different approaches. In Figure 28, the overall 
corneal deformed shapes for the intact and ablated or cut corneal models are compared based 
5 on the fluid-filled porous material assumptions for soft layers. The peripheral portion of the 
cornea alter ablation is thicker than that before ablation. This peripheral thickening matches 
the conceptual model supported by clinical results, described earlier in this application. In 
Figure 29, the overall corneal deformed shapes between the intact and ablated or cut corneal 
models are compared based on the solid material assumptions for soft layers. As opposed to 
10 the results based on the fluid-filled porous material assumption shown in Figure 28, the 
peripheral thickening after ablation is not observed. Therefore, the existing simulations based 
on the solid material assumption are not consistent with clinical results after laser refractive 
surgery, described in this application. 

An ablation depth dependent biomechanical response was described, such that more 
15 ablation implies more flattening. Preliminary calculations are performed for a two-set layer 
ablation where a part of the layer #1,2,3,4 (instead of #1 and #2 as in Figure NK2) is ablated 
or cut. Based on the fluid-filled porous material assumption for soft layers, more flattening is 
observed. However, based on the solid material assumption for soft layers, more bulging is 
observed. This is because the structure rigidity of the central portion of cornea is further 
20 eroded due to further thinning. Therefore, the clinical trend supported in this application is 
reproduced based on the simulations based on the invention, while the wrong trend is 
reproduced by simulations based on the existing methods. 

In addition, current age nomograms for myopic procedures generally tend to require 
less ablation with increasing age, to achieve the same level of correction. It is also known 

25 that the corneas of older individuals tend to be stiffer than those of younger individuals. 
Preliminary parametric studies based on the simple finite element models presented here 
demonstrate that when Young's modulus (layer #1,3,5,7) of the fiber is reduced by a factor of 
10, less flattening in the central portion of the corneal surface is observed. Therefore, the 
trend of basic age nomograms for myopic procedures are reproduced with this innovative 

30 model. 

The importance of this invention is demonstrated through these simple sample 
models. Without the direct incorporation of the negative ("imbibition") pressure inside the 
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cornea into the simulation, the very basic clinically known results relevant to LASIK or PRK 
cannot be properly modeled 

Extension of the sample finite element models 

In some of the existing simulations of refractive surgeries, the swelling due to the 
5 negative "imbibition" pressure in the cornea is indirectly discussed in terms of "hydration" 
(Roy et al. 1996). However, the clear notion that the cornea is highly porous is not directly 
modeled into the formulation. Therefore, there is no separation between the fluid and the 
solid constituents inside the cornea, and the cornea is treated as a ce hydrated" solid material. 
Direct incorporation of negative ("imbibition") pressure and the position (and time)- 
10 dependent change of this pressure due to refractive surgeries into the formulation have never 
been reported in literature. 

This leads to the fact that the key part of this invented modeling can be easily 
extended to all the other existing simulations of refractive surgeries (R. P. Vito, T. J. Shin, 
and B. E.. McCarey, "A Mechanical Model of the Cornea: The Effects of Physiological and 

15 Surgical Factors on Radial Keratotomy Surgery," Refractive & Corneal Surgery, 1989, p. 82- 
88. P. Pinsky and D. V. Datye, "A Microstructrurally-Based Finite Element Model of the 
Incised Human Cornea," Journal of Biomechanics, vol. 24, 1991, p. 907-922. P. Pinsky and 
D V. Datye, "Numerical Modeling of Radial Astigmatic, and Hexagonal Keratotomy, 
Refractive & Corneal Surgery, vol. 8, 1992, p.' 164-172. K. D. Hanna, F. E. Jouve, G. O. 

20 Waring, and P. H. Ciariet, "Computer Simulation of Arcuate Keratotomy for Astigmatism," 
Refractive & Corneal Surgery, vol. 8, 1992, p. 152-163. P. Roy, W. M. Petroll, A.E. 
McKinney an C. J. Chuong, "Computational Models of the Effects of Hydration on Corneal 
Biomechanics and the Results of Radial Keratotomy, Journal of Biomechanical Engineering, 
Transactions of the ASME, vol. 118, 1996, p. 255-258. M. R. Bryant, V. Marchim T. Juhasz, 

25 '^Mathematical Models of Picosecond Laser Keratomileusis for High Myopia," Journal of 
Refractive Surgery, vol, 16, 2000, p, 155-162). As discussed before, in existing simulations, 
the cornea is treated as a homogeneous "hydrated" solid material. All we need to do is to 
make clear throughout the formulation that this homogeneous "hydrated" solid material is in 
fact an "equivalent" homogeneous solid material, and this "equivalent" homogeneous solid 

30 material is actually highly heterogeneous in the sense that it contains pores and pores are 
filled with water. 
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Conceptually speaking, porosity (pore volume per unit volume of a porous solid 
material) can be introduced at each point of this equivalent homogeneous material, and the 
deformation due to the internal pore pressure, which is missing in the existing methods, can 
be considered. Theoretically speaking, the volume expansion due to internal pore pressure 
5 can be simulated by mimicking thermal expansion due to temperature change. Given the 
specific value of pore pressure inside the cornea, this replacement can be done by knowing 
the overall material properties of the porous solid material (solid matrix and pore) and those 
of the solid matrix (without pores). As is explained in the example models, since the cornea 
is highly porous, much of the overall volume change can be attributed to the change in pore 

10 space rather than the expansion or compaction of the solid matrix (without pores). Therefore, 
by employing the simplified assumption of incompressibility of the solid matrix (without 
pores), the expansion due to the given pore pressure can be simulated through thermal 
expansion terms. This leads to the direct incorporation of the negative ("imbibition") pressure 
applied to the pore boundary of the highly porous cornea. In this way, the existing 

15 simulations can be modified to include the effect of the position-dependent negative 
("imbibition") pressure inside the cornea on the corneal biomehanical response to the 
surgery. 

In the simple example models, linearly elastic and isotropic material assumption is 
employed. However, the material response assumption can be extended to anisotropic and/or 

20 nonlinear models. Therefore, by incorporating the basic idea of the existing theories for 
fluid-filled porous materials into the existing simulations, we can significantly improve the 
existing simulations for all other refractive surgeries. 

In the simplified example models, the hard layer representing the collagen fibers and 
the soft layer representing the matrix or ground substance and water are modeled separately. 

25 However, this layered structure can be modeled as an "equivalent" homogeneous anisotropic 
material on the basis of the composite material models for layered structure by Katsube and 
Wu ("A Constitutive Theory for Porous Composite Materials," International Journal of Solids 
and Structures, Vol. 35, pp. 4587-4596, 1998). By introducing the volume fraction of each 
layer (thickness ratio) and the porosity at each point of an "equivalent" homogeneous 

30 material, deformation of the fiber and the ground substance and water inside the cornea can 
be separated within the framework of a single, continuum model. This homogenization 
technique can again be used as a tool to introduce more detailed micro-structure to the 
existing simulations. If the detailed information regarding the microstructure is introduced in 
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the modeling, more detailed information can be obtained from simulations. In addition to the 
above analytical methods, these processes of introducing more microstructure into modeling 
can also be achieved by simply carrying out several levels of macro-micro numerical 
simulations. For example, in the smallest scale level, local positional arrangement 
5 (orientation) of individual collagen fibers may be modeled. In addition, the reorientation and 
stretch of collagen fibers as well as undulation structure of individual collagen fibers can also 
be separately modeled. These small-scale level models can be incorporated into the 
intermediate-scale level of layered structural model. The obtained intermediate-scale level 
model can further be incorporated into the large-scale boundary value problems of corneal 
10 deformation. 

The example models are used simply to demonstrate the key ideas of this invention. 
The same idea can be extended to include various other factors such as time dependent 
movement of fluid and electrochemical considerations. Initial boundary value problems of 
fluid flow can be added to the model simulations so that the time-dependent nature of corneal 

15 deformation can be modeled. In addition to the above purely mechanical consideration of 
cornea, electrochemical balance as well as ion transport can be added to the model 
simulation. The ion and fluid transport phenomena are coupled and both influence the 
deformation of cornea. In both cases, the direct incorporation of negative ("imbibition") term 
becomes very important in solving initial boundary value problems since the atmospheric 

20 pressure that appears at the ablated or cut surface can be specified as boundary conditions. 

The simple finite element models based on this invention can predict the correct trend 
in corneal response to laser refractive surgery, as opposed to existing models in the literature. 
With minor adjustments in geometry and microstructure, the simple models used to describe 
this invention are immediately useful in predicting the general trend in outcomes of refractive 
25 procedures. 

In a LASIK procedure, a flap is first created by a microkeratome and then corneal 
ablation is performed. After ablation, the hinged flap will be replaced in its original position. 
Upon further refinement of the simulations based on this invention, material properties 
relevant to the model can be back calculated from the biomechanical response of a cornea 
30 due to the flap creation (but before ablation). The mechanical response due to ablation can 
then be simulated and the post-operative condition can be predicted based on the model 
simulations. This information can assist the physicians to make decisions about the operation 
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before actual ablation. Therefore, this invention can be very useful in increasing the success 
rate of refractive surgery. 

Integrated System for Refractive Surgery 

The modeling, optimization and algorithm correction of the present invention may be 
5 achieved by written calculation or by implementation via a computer system. Those skilled 
in the art of mathematical computer modeling, computer programming, and database 
manipulation and administration will readily appreciate that the methods and systems 
described herein may be software based and may be performed by a computer system. 

For example, a computer comprising a central processing unit, a storage medium 
10 (such as a magnetic hard drive), a display device (such as a monitor) and an input device 
(such as a keyboard and/or a mouse) may be used to perform the methods of the present 
invention. Pre-operative and post-pertubation measurements are inputted to the computer 
system by keystroke or by other similar methods, including direct data transfer of 
measurements inputted into another computer or similar machine. Existing algorithms for 
15 ablative procedures may be inputted in digital form by any appropriate method. Additional 
empirical, deterministic or other data may be inputted to the computer by any appropriate 
method. 

Data inputted into the computer is stored by any appropriate method, including read- 
only memory (RAM) or other storage devices, such as a hard disk drive. 

20 Analysis and comparison of inputted and/or stored data items is achieved by 

implementing appropriate software algorithms on the computer which compare the relative 
values of stored data. Adjustment of existing algorithms for ablative procedures is achieved 
by retrieving a ablative procedure from data storage and modifying it. The modified 
algorithm is restored or implemented during surgery. 

25 Data derived or created by the present invention is stored in a computer by any 

appropriate method, including storage on an appropriate computer readable medium. An 
example of an appropriate computer readable medium is RAM or a magnetic hard drive. 
Depending upon the nature of the data, the data may be stored in a table or other appropriate 
data structure on such a computer readable medium. Those skilled in the art of databases will 

30 readily appreciate that multiple types of data structures are appropriate under the present 
invention. 
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Modeling is achieved by implementing appropriate software algorithms on the 
computer which use stored data in conjunction with mathematical modeling algorithms. 
Those skilled in the art of computer modeling will readily appreciate that many appropriate 
software modeling applications are available and appropriate under the present invention. 

5 A computer system implementing a method of the present invention may be used 

during cornea surgery. Data may be inputted to such a computer system during the surgical 
process, and the computer system outputs responsive data during the surgery. For example, a 
surgeon or other personnel at surgery inputs pre-operative measurements into the computer. 
The computer, implementing the present invention, determines specifications for ablation 

10 based upon such measurements. After cutting an incision or other pertubation, post- 
pertubation measurements are acquired by diagnostic devices (corneal topography, wavefront 
analyzer, optical coherence tomography, ultrasound) which are integrated with the laser to 
allow all measurements to share a common reference axis. These data are inputted into the 
computer via the integrated system. The computer determines adjusted specifications for 

15 ablation, and outputs such adjustments to the surgeon. The surgeon implements the ablation 
with the adjusted specifications. 

A computer system implementing a method of the present invention may also be used 
prior to cornea surgery to create a biomechanical model based upon inputted data. A 
computer system implementing a method of the present invention may also be used 
20 subsequently to cornea surgery to collect data from the surgery to be used in biomechanical 
modeling. 

What is the Impact of Corneal Biomechanics on Customized Ablative Procedures? 

The studies described above clearly show that corneal biomechanics, including the 
reshaping of the cornea that results from a microkeratome incision, the reshaping of both the 

25 central area of the cornea and the surrounding peripheral area that results from ablative 
surgery, and the reformation produced by post-operative healing must be accounted for to 
establish truly customized ablative procedures for individual patients. Current ablative 
algorithms do not adequately reflect these biomechanical factors. Thus, surgical results to 
date have frequently fallen short of achievable optimum goals. Fortunately, however, 

30 manufacturers of diagnostic and surgical equipment, such as corneal topographers, optical 
coherence tomographers, wave front analyzers and lasers, and surgeons, can account for these 
factors using the procedures and instruments described herein. 20/40 visual acuity with plus 
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or minus one diopter of residual refractive air should no longer be acceptable results for 
LASIK regular or PRK regular surgery. 20/20 or even 20/10 visual acuity, with minimal 
aberration, should be achievable with most patients. 

What is the bottom line for the future of customized laser refractive surgery? Can we 
5 better predict post-operative shape and visual performance? Can we reach toward 20/10 
visual acuity with minimized aberrations? These are absolutely achievable goals if we 
unravel the controllable and/or predictable sources that contribute to the final corneal shape 
and visual outcome. Which instruments will be important to "guide" the laser ablation? It is 
highly likely that both wavefront analysis and corneal topography will be necessary to 

10 program the laser profile in order to consistently achieve an optimized result. Wavefront 
analysis will provide the means necessary to measure and minimize aberrations. Topography 
will help us measure and predict the biomechanical corneal response in ways that have not 
yet been elucidated. The future of customized laser refractive surgery is certainly bright. 
Those skilled in the art of corneal surgery, or the development of diagnostic or surgical 

15 equipment for these procedures will readily appreciate that the methods and instruments 
described herein can produce dramatic improvements in corneal surgery. Of course, they will 
also realize that many modifications may be made in these procedures and instruments, which 
are merely illustrative. These illustrations are not intended to limit the scope of this 
invention, which is defined by the following claims. 
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We claim: 

1 . A method for improving refractive ophthalmic treatment comprising: 
obtaining biomechanical data for the eye from a statistically significant number of 

corneas; and 

establishing a specification for corneal ablation based at least in part on said 
biomechanical properties. 

2. A method in accordance with claim 1 wherein said specification is based at 
least in part on pre-operative data for live corneas. 

3. A method according to claim 1 further comprising perturbing a live cornea to 
generate a biomechanical response in said cornea, and obtaining post-perturbation data 
associated with said biomechanical response. 

4. A method according to claim 3 wherein said perturbation comprises at least 
one of creating a flap on an anterior surface of a cornea, ablating a portion of a cornea, 
peeling the epithelial layer from a cornea, and mechanical or ultrasonic deformation of a 
cornea. 

5. A method in accordance with claim 3 wherein said pre-operative data and said 
post-perturbation data comprises measurement of elevation, total corneal thickness, curvature 
or wave front data. 

6. A method according to claim 5 wherein: 

said perturbation comprises creating a flap on an anterior surface of a cornea; and 
said post-perturbation data comprises flap thickness measurement. 

7. A method in accordance with claim 3 wherein said pre-operative data and said 
post-perrturbation data are obtained by a corneal topographer, an optical coherence 
tomographer, a wave front analyzer, an ultrasonic device, an autorefractor or other diagnostic 
device 
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1 8. A method for performing refractive ophthalmic corneal surgery comprising: 

2 obtaining preoperative diagnostic data that is indicative of at least one biomechanical 

3 property of an eye; 

4 determining ablation specifications based upon said pre-operative diagnostic data; and 

5 ablating a portion of a cornea according to said adjusted ablation specification. 

1 9. A method according to claim 8, further comprising: 

2 perturbing said eye to generate a biomechanical response in said eye; 

3 obtaining post-perturbation diagnostic data of the eye that is indicative of said 

4 biomechanical response; and 

5 determining an adjusted ablation specification based at least in part on a comparison 

6 of said pre-operative data and said post-perturbation data. 

1 10.. A method according to claim 9 wherein said perturbation comprises at least 

2 one of creating a flap on an anterior surface of a cornea, ablating a portion of a cornea, 

3 peeling the epithelial layer from a cornea, and mechanical or ultrasonic deformation of a 

4 cornea. - 

1 1 1 . A method according to claim 9 wherein: 

2 said perturbation comprises creating a flap on an anterior surface of a cornea; and 

3 said post-perturbation data comprise data of flap thickness. 

1 12. A method according to claim 9 further comprising obtaining post-operative 

2 diagnosis data at selected post-operative increments to assess changes in said cornea in the 

3 healing process. 

1 13. A method according to claim 12 wherein said post operative increments 

2 comprise at least one of about one day, about one week, about one month, about six months, 

3 and about nine months following said ablation. 

1 14. A method in accordance with claim 12 comprising phdtoablative retreatment 

2 of said cornea. 
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15. A method for establishing a specification for photoablative ophthalmic surgery 
comprising: 

obtaining pre-operative ophthalmic data that is indicative of at least one 
biomechanical response of a live cornea; and 

determining ablation specifications based upon said pre-operative data. 

16. A method according to claim 1 5 further comprising: 
perturbing a cornea to generate a response; 

obtaining post-perturbation ophthalmic data; and 

determining an adjusted ablation specification based at least in part on a comparison 
of said pre-operative diagnostic data and said post- perturbation data wherein these 
measurement are indicative of a biomechanical response . 

17. A method for improving refractive ophthalmic surgery comprising: 
obtaining pre-operative data for a statistically significant number of live corneas; 
predicting biomechanical response to ablation of said corneas; and 

establishing one or more mathematical models of the biomechanical response of said 
corneas to said perturbations. 

18. A method according to claim 17 wherein said mathematical models are based 
on finite element analysis or other mathematical modeling techniques. 

19. A method according to claim 17, further comprising: 

perturbing a statistically significant number of said corneas to generate biomechanical 
responses in said corneas; 

obtaining post-perturbation diagnostic data for said corneas that is indicative of said 
biomechanical response; and 

establishing said mathematical models based at least in part on a comparison of said 
pre-operative data and said post-perturbation data. 

20. . A method according to claim 1 9 wherein perturbation comprises at least one 
of creating a flap on an anterior surface of a cornea, ablating a portion of a cornea, peeling the 
epithelial layer from a cornea, and mechanical or ultrasonic deformation of a cornea. 
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2 2 1 . A method according to claim 1 7 further comprising: 

3 ablating a portion of said statistically significant number of corneas; 

4 obtaining post-operative diagnostic data for said statistically significant number of 

5 said corneas: and 

6 determining an adjusted ablation specification based at least in part on a comparison 

7 of said pre-operative data, said post- perturbation data and said post-operative data, wherein 

8 said comparison is indicative of a biomechanical response . 

1 22. A method according to claim 21 comprising refractive retreatment of one or 

2 more of said corneas. 

1 23. A method according to claim 17 wherein the process for establishing said 

2 mathematical model comprises presenting said model as a portion of a hydrated porous 

3 structure. 

1 24. A method according to claim 23 wherein said structure comprises a plurality 

2 of hard layers separated by at least one soft, porous layer; and collagen fibers are modeled as 

3 impermeable hard layers. 

1 25. A method according to claim 24 wherein said collagen fibers are collectively 

2 assumed to be thin hard shell layers. 

1 26 A method according to claim 23 wherein said cornea is modeled as a layered 

2 or lamelar, highly porous structure, comprising a plurality of said thin hard shell layers 

3 separated by at least one soft porous layer. 

1 27. A method according to claim 26 wherein said porous material comprises a 

2 matric of incompressible substance with pores inside said matrix, at leawst some of said 

3 pores being at least partially filled with fluid. 

1 28. A method for predicting biomechanical response of live corneas, comprising: 

2 obtaining data from a statistically significant number of corneas; 

3 creating a finite element model of a cornea, comprising: 
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4 a portion of a hydrated porous shell: 

5 fixed edges; 

6 a substantially uniform intraocular pressure of about 10 to 

7 about 15 mm of Hg against an inner surface of said shell: and 

8 a subatmospheric pressure within said shell; and 

9 simulating a localized change in the fluid content of said shell; and 

1 0 measuring a biochemical response of said shell following said change. 

1 29. A method for correcting an existing algorithm for ablative procedures in 

2 cornea surgery comprising: 

3 receiving pre-operative data that is indicative of at least one biomechanical response . 

4 of a live cornea; and 

5 adjusting said existing algorithm for ablative procedures based upon said pre- 

6 operative biomechanical data. 

1 30. A method according to claim 29 further comprising: 

2 perturbing said cornea: 

3 receiving post-perturbation data of the cornea before ablating said cornea 

4 analyzing said pre-operative data and said post-perturbational data; and 

5 adjusting said existing algorithm for ablative procedures based upon differences 

6 between said pre-operative data and said post-perturbation data. 

1 3 1 . A method according to claim 30, further comprising: 

2 ablating said cornea in accordance with said adjusted algorithm; 

3 taking post-operative data of said cornea; and 

4 establishing corrections to said existing algorithms for future surgery based on 

5 differences between said differences between said pre-operative data and said post-operative 

6 data. 

1 32. A method for establishing a biomechanical model for customized ablative 

2 procedures in refractive cornea surgery for individual patients comprising: 

3 receiving pre-operative data for a statistically significant number of live corneas; 

4 perturbing said corneas to generate a biomechanical response; 

5 receiving post-perturbation data for said corneas; and 
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6 establishing one or more models of the biomechanical response, of said corneas to said 

7 perturbations. 

1 33. A method according to claim 32 further comprising ablating anterior surfaces 

2 of said corneas, receiving post-operative data for said corneas following said ablation, and 

3 establishing one or more models for the biomechanical reponse of said corneas to said 

4 perturbations and said ablation. 

1 34. A method according to claim 32 further comprising receiving post-operative 

2 healing data, and establishing one of more models for the biomechanical response of said 

3 corneas to said perturbations and said healing. 

1 35. A method according to claim 32 wherein said data is received in a computer 

2 system, which generates said model or models, for establishing specifications for ablative 

3 procedures in refractive cornea surgery comprising: 

1 36 A computer-readable medium containing a data structure for storing 

2 corrections for existing algorithms for ablative procedures in cornea surgery, comprising: 

3 an algorithm table containing at least one entry for each of a plurality of existing 

4 algorithms for ablative procedures, 

5 each said entry containing a plurality of corrections to said existing algorithm, 

6 each said plurality of corrections containing corrective specification for said existing 

7 algorithm based upon pre-operative data of a cornea and post-perturbational data for a cornea. 

1 3 7. An integrated system for performing refractive cornea surgery comprising: 

2 means for supporting a patient in a therapeuticposition and thereby establish a 

3 common reference axis via fixation point with respect to a cornea of said patient; 

4 a surgical instrument to provide a flap on an anterior surface of a cornea and means to 

5 position said surgical insrument with respect to said axis; 

6 one or more measuring instruments adapted to be positioned with respect to said 

7 reference axis for taking data of said cornea, while said patient is supported in said 

8 therapeutic position, before and after said flap has been formed; and 

9 a laser adapted to be positioned with respect to said reference axis instrument for 

10 ablating said cornea while said patient is in said supported in said therapeutic position.. 
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1 38. An integrated system in accordance with claim 37 wherein the instrument for 

2 providing said flap comprises a keratome or a laser. 

1 39. An integrated system in accordance with claim 37 said pre-operative data and 

2 said post-perturbation data are obtained by a corneal topographer, an optical coherence 

3 tomographer, a wave front analyzer, an ultrasonic device, an autorefractor or other diagnostic 

4 device 
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Figure 1: Schematic of the simple "shape^subtraction" 
paradigm for correction of myopia. K\ and Ra are initial ^ 
are final radii of curvature; to is the maximum depth of the 
laser cut; S is the diameter of the optical zone. Reprinted 
from Munnerlyn et al. 




Figure 2: Schematic diagram of shape-subtraction model of refractive surgery for a myopic 
ablation (left) and hyperopic ablation (right). The pre-operative radius of curvature is Rl and the 
desired post-operative curvature is R2. The mayi'miim ablation depth is det ermine d by the ablation 
zone size. The intervening tissue, shown in red, between the pre-operative curve (solid) and post- 
operative curve (dashed/dotted) is "subtracted" with an excimer laser to produce the desired result. 
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Figure 3: Post-operative minus pre-operative elevation (upper left), pachymetric (lower left), 
tangential curvature (upper right) and axial (lower right) difference maps, measured with an 
Orbscan II for a patient with a 6mm ablation zone. The white circles define three corneal 
regions. The central region has radius 2.75mm, the middle region has radius 2.75 to 3.25mm and 
the outer region has radius 3.25 to 4.5mm. Elevation (upper left) differences were calculated by 
fitting the two surfaces in the middle region. Note the increases in elevation, pachymetry and 
tangential curvature in the region outside the ablation zone. 
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Figure 4: Tahgetitial clirvaturd difference map (post-op minus pre-op) aftet a -12.5 diopter 
LAS DC procedure With a 5.5rnm diameter ablation zone. White pifcels indicate zero difference. 
Central red circle is at 4mm diameter and outer white circle is at 8mm diameter. Note the 
increased curvature outside the ablation zone. 
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Figure 5: Pachymetry difference map (post-op minus pre-op) after a -12.5 diopter myopic 
LASIK procedure with a 5.5mm diameter ablation zone. White pixels indicate zero difference. 
Central red circle is at 4mm diameter and outer white circle is at 8mm diameter. Note the 
increased pachymetry outside the ablation zone. 
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Figure 6: Post-op minus pre-op difference maps from a -2.25 diopter LASIK correction with a 
6mm diameter ablation zone. Topography was acquired using Orbscan and exported for analysis 
using custom software. The inner, middle and outer white circles have radii of 2.75mm, 3.25mm 
and 4.5mm, respectively. The central zone within the inner circle lies within the ablation zone. 
The edge of the ablation zone is contained within the annulus between the center and middle circle. 
Expected central decreases in elevation, pachymetry, and curvature are seen as negative 
differences on the respective maps. However, unexpected increases in elevation, pachymetry and 
curvature have occurred outside the ablation zone, due to the biomechanical response of the cornea 
to a change in its. structure. This is seen within the annulus between the middle and outer circles, 
where the redder colors represent increases. 
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Figure 7: Error maps of actual minus predicted post-op topography for the same subject as Figure 
6. Predicted post-op topography was generated by calculating an estimated ablation profile using 
Munuerlyn's formulas and subtracting it from the pre-operative raw elevation. Predicted curvature 
was then calculated along corneal meridians. Note that both the elevation and pachymetry error 
maps demonstrate less than predicted central values (or greater than expected ablation depth), as 
indicated by negative numbers and blue color. However, both maps also paradoxically show 
greater than predicted differences outside the ablation zone. The tangential curvature error map 
also shows greater than expected curvature change outside the ablation zone. Centrally, the values 
range from less than predicted to greater than predicted, with overall less error than outside the 
ablation zone. These errors are consistent with the biomechanical model proposed herein. 
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Figure 8: A) Pre-operative elevation (left) and tangential curvature (right). B) Post-onerative 

ir S } aDd ?lf tia ! curvature (right) -*«■ a - 7 - 63 pfodST % 

mevation error map Oeft) and tangential curvature error map (right) generated by subtotini 
known ablation algorithms (supplied by Summit TechnologS) f LTufe L oSvf 
topo^phy (A) and generating a predicted post-operative topography, which w^Eed 
from the measured topography (B) to create the difference maps X red S 

Sa C ™S^ ential , eiror r P repreSOTtS 1,16 central P-treamlent anTfoe ou£ red 
area represents the biomechanical response of the cornea. 
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Figure 9: A) Pre-operative elevation (left) and tangential curvature (right). B) Post-operative 
elevation (left) and tangential curvature (right) after a -3.82 diopter LASIK procedure. C) 
Elevation error map (left) and tangential curvature error map (right) generated by subtracting 
known ablation algorithms (supplied by Summit Technologies) from the pre-operative 
topography (A) and generating a predicted post-operative topography, which was subtracted 
from the measured topography (B) to create the difference maps shown in (C). The red area in 
the center of the tangential error map represents the central island pre-treatment and the outer red 
area represents the biomechanical response of the cornea. 
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Figure 10: Regression analysis of peripheral stromal thickness of the superior region (upper plot) 
and inferior region (lower plot) against curvature in the central 3 mm region, pre-operatively and at 
4 time points postoperatively after PRK in one patient. This demonstrates how central curvature 
for this patient closely tracked peripheral stromal thickness over time. 
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Figure 11: A conceptual model that predicts biomechanical central flattening as a direct 
consequence of severed corneal lamellae. (A- TOP) The pre-operative cornea may be conceived 
as a senes of stacked rubber bands (lamellae) with sponges between each layer (interlamellar 
spaces rilled with ground substance or matrix). The rubber bands are in tension, since there is a 
force pushing on them from underneath (intraocular pressure), and the ends are held tightly by 
SS^SfT" ^ ermy ° pic lase "«***ve surgery, a series of lamellae areTeverS 

and obliterated centrally The remaining peripheral segments relax, just like the taut rubber 
bands would relax once cut. With the reduction of tension in the peripheral lamellae, the 
squeezing force on the matrix is reduced and the distance between lamellae expands, analogous 
to the sponges taking up water if the rubber bands are cut This allows the periphery of the 
cornea to thicken. Due to the crosslinking between lamellar layers, the expansion force nulls on 
die underlying intact lamellae, as indicated by the arrows pointing radially outward. An outward 
force in the periphery pulls laterally on the center and flattens it. Thus, the cornea will flatten 
centrally with any procedure that circumferentially severs lamellae. This includes myopic 
profits, hyperopic profiles, constant depth-PTK profiles, as well as the simple cutting of a 
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Interlamellar COHESIVE Strength 



& 



Interlamellar SHEAR Strength 



Figure 12: A schematic diagrarn of the forces described by A) the interlarnellar cohesive strength 
and fi) the ipterlawell^r shear strer+gtb. 
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BIOMECHANICS EFFECT 
flattening 
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Figure 13: Proposed impact of biomechanics on LASIK nomograms for primary and secondary 
hyperopia. Biomechanical flattening opposes corneal steepening in attempted corrections of 
hyperopia. This effect iS more pronounced in previously untreated corneas (1° hyperopia, A) and 
requires compensatory enhancement of the ablation profile to achieve the desired steepening 
effect. Corneas of eyes that are hyperopic after previous surgical overcorrection (B) are less 
"loaded" biomechanically; i.e., peripheral expansion (blue x's) and associated central flattening 
have already occurred. Accordingly, profile effects are opposed to a lesser degree and little or no 
enhancement of the ablation profile is required. Corneal geometry is exaggerated for illustrative 
purposes. 
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6 Months 





figure 14: Sample ablation depth profile calculated at 6months after PRK by subtracting pre-oo 
2?«I7^f ■ ^ ^ ^ 0f&ttt mCaSUred * a PAR Comea I Topography System. Maximal 
- S ^ * ~ C6nter ° fthe P rofile near * e ^ The edges of the ablation zone are 
hnSl^f^!-. 111 , Pf ^f ? "W"™**"^ ±3mm from the apex at the intersection with the 
Si. ^ ^ d6pictS i relative 2610 0uteide Elation zone, the depth 

t ,1 rel3tlVe 2ero ^ a "^Sative" depth of ablation or mcreased 

2S3i£^ d ^ ° f ^ —Ptptay is consistent with the proposed 
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Figure IS: Comparison between planned ablation profiles using CEPTA (see Chant* yy% «« 
LEFT and measured topographic elevation difference ma^onfe ]ff T ^ 

nS^"° P ^u Pn5 "°? SUrfkces Where ^ ™» rcm °ved, and yellow/ri Shis S 
d^rence where tissue was "added" after surgery. Zero difference HSSrfif ™ 

l™£r£r? l V n 5? kft 216 SUn0Unded ^ ^ zon « indicating zero^™ o^oE' 
zZT^- SV!ltl<m ^r aD0 Sh ° W P ° Sitive e,evation of thTcSrafSon 
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^ r ! U L Av ^ ge Differe,lce Maps between repeated measures of 20 eyes of 10 normal 

£Sn ^ CkCl ^ d6fine C ° meal re *^ ^ ^tral region 2TnSL 2 Sn 
foe middle region has radius 2.75 to 3.25mm and the outer region bu indto 325 to 4 S 

^ASSti^ff^^ diff — dower lenfand 
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Figure 17: Average difference maps between 1 day post-op LASIK and pre-op state for 8 eyes of 
4 patients. The white circles define three corneal regions. The central ablated region has radius 
2.75mm, the middle transitional region has radius 2.75 to 3.25mm and the outer non-ablated region 
has radius 3.25 to 4.5mm. Elevation (upper left) differences were calculated by fitting the two 
surfaces in the transitional region, and demonstrate increased elevation in the outer zone. Average 
tangential differences (upper rigfct) show increased curvature in the outer region. Average 
pachymetry differences (lower left) show increased pachymetry in the outer region after surgery. 
Average axial differences (lower right) are also given. 
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Figure 18: Average difference maps between 1 month post-op LASK and pre-op state for 8 
eyes of 4 patients. The white circles define three corneal regions. The central ablated region has 
radius 2.75mm, the middle transitional region has radius 2.75 to 3.25mm and the outer non- 
ablated region has radius 3.25 to 4.5mm. Elevation (upper left) differences were calculated by 
fitting the two surfaces in the transitional region, and demonstrate persistent increased elevation 
in the outer zone. Average tangential differences (upper right) show persistent increased 
curvature in the outer region. Average pachymetry differences (lower left) show persistent 
increased pachymetry m the outer region. Average axial differences (lower right) are also given 
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Figure 20. Examples of corneal optical section images. The central region (A) and inferior (B) 
and superior (C) peripheral regions constitute a section of the vertical meridian in an ablated 
cornea. Ink dots and corneal landmarks facilitated analysis of corresponding regions in paired 
controls and replicate images. Digital tracings of the anterior and posterior stromal surfaces (not 
shown) were used in conjunction with mathematically-derived projections from the posterior 
corneal surface (shown only in A) to obtain apparent thickness measurements. Corrections for 
refraction and viewing angle were then applied. Central cornea images (A) were subdivided for 
analysis of differential thickness losses across the ablation zone. 
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A. Near Inferior B. Far Inferior 




Thickness Change (urn) Thickness Change (urn) 



C. Near Superior D. Far Superior 




Thickness Change (urn) Thickness Change (urn) 



Figure 21. Anterior curvature shifts were correlated to stromal thickness changes in peripheral 
subregions of the vertical meridian. Negative curvature changes indicate corneal flattening. 
Correlation coefficients were significant in all but the far inferior periphery (P values = .03, .05, 
.03 and .01 for A to D, respectively). 
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Figure 22: Pre-operative Keratron Scout curvature topography taken while the patient is supine 
for a LASIK procedure is given on the upper left. The post-flap, pre-ablation topography is 
given on the lower left This was acquired immediately after cutting the flap, and before it was 
reflected. Note that the characteristic post-refractive surgery pattern of central flattening and 
steep mid-peripheral "red ring" begins to appear after the flap and before the ablation. Finally, 
the post-ablation topography is given on the upper right, acquired after the flap is replaced, and 
with the patient supine. 
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Figure 23: On the left is the curvature difference map between pre-op and post-flap for the 
patient in Figure 25, calculated and displayed using a customized software. The average dioptric 
difference in the central region is -1.08 diopters (flattening), and the average dioptric difference 
in the outer region is +10.20 diopters (steepening), all induced just by cutting the flap. On the 
right is the curvature difference map between post-flap and post-ablation for the patient in Figure 
5. The average dioptric difference in the central region is -2.25 diopters (further flattening), and 
the average dioptric difference in the outer region is +8.37 diopters (further steepening). 
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Figure 24 

Intact corneal model subjected to the intraocular pressure (15 mmHg) at the 
posterior surface and the negative pressure ("Imbibition pressure", 60 mmHg) at 
the pore boundary of the highly porous soft layers #2, 4 and 6. 



Axisymrnetric 




A schematic presentation of a highly porous material 
subjected to the negative pressure ( "imbibition pressure", 60 mmHg ) 
applied to the pore boundary 
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Figure 25 

Cut corneal model subjected to the intraocular pressure (15 mmHg) at the posterior 
surface and the negative pressure ("Imbibition pressure", 60 mmHg) at the pore 
boundary of the highly porous soft layers #4 and 6; zero pore pressure for soft 
layer #2. 
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Figure 28 

Comparison of the overall corneal deformed shape between the intact and cut 
corneal models With highly porous material assumption in soft layers #2,4 and 6: 
Pink region represents the intact cornea 
Blue region represents the cut cornea 
(Displacement scale x10) 
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Figure 29 

Comparison of the overall corneal deformed shape between the intact and cut 
corneal models with solid material assumption (without pores) in soft layers #2,4 
and 6: 

Pink region represents the intact cornea 
Blue region represents the cut cornea 
(Displacement scale x80) 
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